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ABSTRACT 



We present the deepest 100 to 500 yum far-infrared observations obtained with the Herschel Space Observatoiy as part of the GOODS-Herschel 
key program, and examine the infrared (IR) 3-500 yum spectral energy distributions (SEDs) of galaxies at < z < 2.5, supplemented by a local 
reference sample from IRAS, ISO, Spitzer and AKARI data. We determine the projected star formation densities of local galaxies from their radio 
and mid-IR continuum sizes. 

We find that the ratio of total IR luminosity to rest-frame 8 fim luminosity, IRS (=L}^/Lg), follows a Gaussian distribution centered on 7^8=4 
(cr= 1.6) and defines an IR main sequence for star-forming galaxies independent of redshift and luminosity. Outliers from this main sequence 
produce a tail skewed toward higher values of IRS. This minority population (<20 %) is shown to consist of starbursts with compact projected 
star formation densities. IRS can be used to separate galaxies with normal and extended modes of star formation from compact starbursts with 
^ high-/7?8, high projected IR surface brightness (Sir>3x 10'" Lakpc"^) and a high specific star formation rate (i.e., starbursts). The rest-frame, UV- 
5-H 2700 A size of these distant starbursts is typically half that of main sequence galaxies, supporting the correlation between star formation density 
jy^ and starburst activity that is measured for the local sample. 

Ci Locally, luminous and ultraluminous IR galaxies, (U)LIRGs (Lj^>10"Lo), are systematically in the starburst mode, whereas most distant 
(U)LIRGs form stars in the "normal" main sequence mode. This confusion between two modes of star formation is the cause of the so-called 
"mid-IR excess" population of galaxies found at z>1.5 by previous studies. Main sequence galaxies have strong polycyclic aromatic hydrocarbon 
(PAH) emission line features, a broad far-IR bump resulting from a combination of dust temperatures (rdi,st~15 - 50 K), and an effective Tjast~a 
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' K, as derived from the peak wavelength of their infrared SED. Galaxies in the starburst regime instead exhibit weak PAH equivalent widths and 
a sharper far-IR bump with an effective rdust~40 K. Finally, we present evidence that the mid-to-far IR emission of X-ray active galactic nuclei 
. (AGN) is predominantly produced by star formation and that candidate dusty AGNs with a power-law emission in the mid-IR systematically occur 
in compact, dusty starbursts. After correcting for the effect of starbursts on IRS, we identify new candidates for extremely obscured AGNs. 

Key words. Galaxies: evolution - Galaxies: active - Galaxies: starburst - Infrared: galaxies 
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T— I ' 1. Introduction larger at z~l, i.e., ~ 8 Gyr ago, when LIRGs were responsible 

L',-. ,-, , for most of the cosmic SFR density per unit co-moving volume 

_ > , It IS now well established that -85 % of the baryon mass ^^^^ ^ ^1^^^ 2001 - hereafter CEOl, Le Floch et al. 

tX! ■ 'r.'; ^^/'T u ""y^^'^Ii 2005, Magnelli et al. 2009). Eariier in the past, at z~2, sub-mm 

^ ■ Marchesmi et al. 2009 and references therem) and that most en- ^ .^^^^ observations revealed that the contribution to the cos- 

ergy radiated dunng this epoch by newly formed stars was heav- ^.^ ^^^^ ^^^^ ^^^^^ ^^.^^^^^ ultraluminous 

- ■ ■ ily obscured by dust. To understand how present-day galaxies j^. j^^.^^ (ULIRGs, L,^>10'^ Lq), was as important as for 

were niade, it is therefore imperative to accurately determine ^irGs (Chapman et al. 2005, Papovich et al. 2007, Caputi et al. 

the bolometric output of dust hence the total IR luminosity, 2007, Daddi et al. 2007a, Magnelli et al. 2009, 201 1). However, 

L5°', integrated from 8 to 1000 ^m. In the past, this key infor- ^^^^ ^^^^^ ^^^^.^^ ^^^^ rest-frame far-IR measurements of 

mation on the actual star formation rate (SFR) experienced by individual galaxies at wavelengths where the IR spectral energy 

distant galaxies was determined by extrapolating obsei-vations distribution (SED) of star-forming galaxies is known to peak, 

in the mid-IR and sub-milhmeter (sub-mm) or by correcting ^^^^^ ^.^j^^d on stacking of fai--IR data from individually 

their UV luminosities for extinction. These extrapolations im- undetected sources 

plied that the number density per unit comoving volume of lu- , , , ,o ^, /t,„ 

minous IR galaxies (LIRGs, 1O"<Lir/L0<1O12) was 70 times ^ith the launch of the Herschel Space Observatoiy (Pilbratt 

et al. 2010), it has now become possible to measure the total IR 

* Herschel is an ESA space observatory with science instruments luminosity of distant galaxies direcdy. Using shallower i/e«c/!e/ 

provided by European-led Principal Investigator consortia and with im- data than the present study, Elbaz et al. (2010) showed that ex- 

portant participation from NASA. trapolations of from the mid-IR (24 fim passband), which 
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was done under the assumption that the IR SEDs of star-forming 
galaxies remained the same at all epochs, were correct below 
z<1.3, with an uncertainty of only 0.15 dex. However, the ex- 
tension of this assumption to (U)LIRGs at z>1.3, in large part 
relying on stacking, failed by a factor 3-5 typically (Elbaz et 
al. 2010, Nordon et al. 2010). This finding confirmed the past 
discovery of a so-called "mid-IR excess" population of galaxies 
(Daddi et al. 2007a, Papovich et al. 2007, MagnelH et al. 2011): 
the 8 i-im rest-frame emission of z~2 (U)LIRGs was excessively 
strong compared to the IR SED of local galaxies with equivalent 
luminosities when deriving Lj°' from the radio continuum at 1 .4 
GHz, from stacked measurements from Spitze r-MlPS lOfim, or 
from the UV luminosity corrected for extinction. 

Various causes have been invoked to explain this "mid-IR 
excess" population: (i) an evolution of the IR SEDs of galax- 
ies; (ii) the presence of an active galactic nucleus (AGN) heat- 
ing dust to temperatures of a few 100 K; or (Hi) limitations in 
local libraries of template SEDs, i.e., the A:-correction effect on 
distant galaxies probing regimes where the SEDs were not ac- 
curately calibrated. Evidence pointing toward an important role 
played by obscured AGN to explain these discrepancies (point 
(//)) came from the stacking of Chandra X-ray images at the 
positions of the most luminous z~2 BzK galaxies (Daddi et 
al. 2007a). The most luminous of these distant galaxies were 
detected in both the soft (0.5-2 keV) and hard (2-7 keV) X- 
ray channels of Chandra and exhibited a flux ratio typical of 
heavily obscured {Nh>10^^ cm"^) or even Compton thick AGN 
(Nh>^0^* cm"^). Surprisingly, however, a high fraction of the 
same objects, when observed in mid-IR spectroscopy with the 
Spitzer IR spectrograph (IRS), were found to possess intense 
polycyclic aromatic hydrocarbon (PAH, Leger & Puget 1984, 
Puget & Leger 1989, Allamandola et al. 1989) broad lines with 
equivalent widths strongly dominating over the hot to warm dust 
continuum (Rigby et al. 2008, Farrah et al. 2008, Murphy et al. 
2009, Fadda et al. 2010, Takagi et al. 2010). Deeper Chandra ob- 
servations have since showed that only ~25 % of the z~2 BzK- 
selected mid-IR galaxies hosted heavily obscured AGN, the rest 
being otherwise composed of relatively unobscured AGNs and 
star-forming galaxies (Alexander et al. 2011). This would in- 
stead favor points (/) or above. 

In this paper, we present the deepest 100 to 500 yum far-IR 
observations obtained with the Herschel Space Observatory as 
part of the GOODS-Herschel Open Time Key Program with the 
PACS (Poglitsch et al. 2010) and SPIRE (Griffin et al. 2010) in- 
struments. Thanks to the unique power of Herschel to determine 
the bolometric output of star-forming galaxies, we demonstrate 
that incorrect extrapolations of Lj°' from 24 yum observations at 
z>1.5, and the associated claim for a "mid-IR excess" popula- 
tion, do not indicate a drastic evolution of infrared SEDs, nor 
the ubiquity of warm AGN-heated dust dominating the mid-IR 
emission. Instead, we show that the 8 ^m. bolometric correction 
factor {IR%=Lf^lL^) is universal in the range 0<z<2.5, hence 
defining an IR "main sequence" (MS). We show that past incor- 
rect extrapolations resulted from the confusion between galaxies 
with extended star formation and those with compact starbursts, 
which exhibit notably different infrared SEDs. 

We present evidence that this IR main sequence is directly 
related to the redshift dependent SFR - M* relation (Noeske et 
al. 2007, Elbaz et al. 2007, Daddi et al. 2007a, 2009, Pannella 
et al. 2009, Magdis et al. 2010a, Gonzalez et al. 2011) and is 
able to separate galaxies between those experiencing a "normal" 
mode of extended star formation and starbursts with compact 
projected star formation densities. This distinction between a 
majority of "main sequence" (MS) galaxies and a minority of 



compact "starbursts" (SB) is analogous to the recent finding of 
two regimes of star formation in the Schmidt-Kennicutt (SK) 
law, with MS galaxies following the classical SK relation while 
the SFR of SB galaxies is an order of magnitude greater than 
expected from their projected gas surface density (Daddi et al. 
2010, Genzel et al. 2010). To separate these two star-formation 
modes, the GOODS-//eric/ie/ observations of distant galaxies 
are supplemented by a reference sample of local galaxies using 
a compilation of data from IRAS, AKARI, Spitzer, SDSS and ra- 
dio observations. 

The GOOT)S>-Herschel observations and catalogs are pre- 
sented in Section |2] The main limitation of the Herschel cata- 
logs, the confusion limit, and a "clean index" identifying sources 
with robust photometry are discussed in Sect. 12.31 The high- and 
low-redshift galaxy samples are introduced in Section|3]together 
with a description of the method used to compute total IR lumi- 
nosities, stellar masses and photometric redshifts. The IR main 
sequence is presented in Section |4] where the so-called "mid-IR 
excess problem" is addressed and a solution proposed using the 
IR% bolometric correction factor This parameter, which relies on 
the same rest-frame wavelengths independent of galaxy redshift, 
is used to separate star-forming galaxies in two modes: a main 
sequence and a starburst mode. In the following sections, IR% is 
shown to correlate closely with the IR surface brightness, hence 
with the projected star formation density, and with the starburst 
intensity, that we quantify here with a parameter named "star- 
burstiness", for local (Section |5]l and distant (Section |6]l galax- 
ies. It is shown that galaxies exhibiting enhanced IR% values are 
undergoing a compact starburst phase. The universality of //?8 
among main sequence star-forming galaxies is used to produce 
a prototypical IR SED for galaxies in the main sequence mode 
of star formation in Section|7] We combine Spitzer and Herschel 
photometry in many passbands for galaxies at < z < 2.5 to 
derive composite SEDs for both main sequence and starburst 
galaxies. Finally, galaxies exhibiting an AGN signature are dis- 
cussed in Section [8] where we present a technique to identify 
obscured AGN candidates that would be unrecognized by previ- 
ous methods. 

We use below a cosmology with //o=70 kms^'Mpc"', Q.m - 
0.3, Qa = 0.7 and we assume a Salpeter initial mass function 
(IMF, Salpeter 1955) when deriving SFRs and stellar masses. 

2. GOODS-Herschel data and catalogs 

2.1. Observations 

The sample of high-redshift galaxies analyzed here consists 
of galaxies observed in the two Great Observatories Origins 
Deep Survey (GOODS) fields in the Northern and Southern 
hemispheres. Observations with the Herschel Space Observatory 
were obtained as part of the open time key program GOODS- 
Herschel (PI D.Elbaz), for a total time of 361.3 hours. PACS 
observations at 100 and 160 yum cover the whole GOODS-north 
field of 10'xl6' and part of GOODS-south, i.e., lO'xlO' (but 
reaching the largest depths over ~64 arcmin^). When consider- 
ing the total observing times of 124 hours in GOODS-N and 
206.3 hours in GOODS-S (including 2.6 and 5 hours of over- 
heads), the PACS GOODS-Herschel observations reach a total 
integration time per sky position of 2.4 hours in GOODS-N and 
of 15.1 hours in GOODS-S, i.e., 6.3 times longer Due to the 
larger beam size and observing configuration, the SPIRE obser- 
vations of GOODS-N cover a field of 900 arcmin^, hence largely 
encompassing the central 10'xl6', for a total observing time of 
31.1 hours and an integration time per sky position of 16.8 hours. 
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Fig. 1. Postage stamp images of the same 5'x5' region of the GOODS-north field ranging from 3.6 yum (upper-left) to 500 jum 
(bottom-right). Upper panel: five Spitzer images (85cm telescope diameter) obtained with IRAC at 3.6 and 8 /vm, the IRS peak-up 
array at 16/^m and MIPS at 24 and 70/im (from upper-left to upper-right). Bottom panel: five Herschel images (3.5m telescope 
diameter) obtained with PACS at 100 and 160 yum and SPIRE at 250, 350 and 500 //m (from bottom-left to bottom-right). 




Fig. 2. Composite three color image of the GOODS-north field 
(10'xl5') at 100 A/m (blue), 160 jum (green) and 250jum (red). 
North is up and east is left. 



Fig. 3. Composite three color image of the GOODS-south field 
(lO'xlO') at 24fim (blue), 100 jum (green) and 160 //m (red). 
North is up and east is left. 



Fig.[T]shows a montage of images (each 5'x5') from Spitzer- 
IRAC at 3.6yum to SPIRE at 500 yum. This illustrates the impact 
of the increasing beam size as a function of wavelength: the 
number of sources that are clearly visible at each wavelength 
increases when going from the longest to the shortest wave- 
lengths (with the exception of the 70;um image, which comes 
from Spitzer and not Herschel). Composite three color images 
of GOODS-N at 100-160-250yum and GOODS-S at 24-100- 
160 yum are shown in Figs.|2]and[3] 



2.2. Catalogs 

2.2.1. Source extraction 

Flux densities and their associated uncertainties were obtained 
from point source fitting using 24 /im prior positions. For the 
largest passbands of SPIRE (i.e., 350 and 500 //m), the 24 //m 
priors are much too numerous and would lead to an over- 
deblending of the actual sources. Hence, we defined priors with 
the following procedure. For PACS-lOOyum, we used MIPS- 
24 yum priors down to the 3cr limit and imposing a minimum flux 
density of 20//Jy. For PACS-160jum and SPIRE-250yum, we re- 
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Fig. 4. GOODS-Herschel detection limits (from Table [TJ 
and total IR luminosities of the Herschel sources as a 
function of redshift (filled dots: spectroscopic, open dots: 
photometric). The right axis is the SFR derived from 
SFR[Moyr-i] = 1.72xlO-i"xL™[L0] (Kennicutt 1998a). These 
limits were computed assuming the local library of template 
SEDs of CEOl. For comparison, the Spitzer MIPS-24/vm de- 
tection limit is represented as well as the knee of the total IR 
luminosity function, as derived by Magnelli et al. (2009, 2010). 

striated the 24 /im priors to the 5cr (SOyuJy) hmit (reducing the 
number of priors by about 35 %). For SPIRE-350 and 500 ;um, 
we kept only the 24 yum priors for sources with a S/N ratio greater 
than 2 at 250 yum. These criteria were chosen from Monte Carlo 
simulations (see Sect. 12.2.21 ) to avoid using too many priors that 
would result in subdividing flux densities artificially, while pro- 
ducing residual maps (after PSF subtracting the sources brighter 
than the detection limit) with no obvious sources remaining. 

Fig. m can be used to infer the reliability of the Spitzer- 
MIPS 24nm images of the two GOODS fields for identifying 
potential blending issues with Herschel. It shows that the 20 juJy 
depth at 24 /vm (3cr) reaches fainter sources than any of the 
Herschel bands, down to the confusion level and up to a red- 
shift of z~3. The technique used to estimate total IR luminosi- 
ties for the Herschel sources is discussed in Sect. 13.21 Hence the 
positions of 24 /im sources can be used to perform robust PSF 
fitting source detection and flux measurements on the Herschel 
maps. We validated the efficiency of this technique by checking 
that no sources remain in the residual images after subtracting 
the detected sources or by independently extracting sources us- 
ing a blind source extraction technique (Starfinder: Diolaiti et al. 
2000). Although it is the case that most Herschel sources have 
a 24 jum counterpart, a few 24 //m-dropout galaxies were found, 
i.e., galaxies detected by Herschel but not at 24yum. This will be 
the subject of a companion paper (Magdis et al. 201 1). But these 
objects represent less than 1 % of the Herschel sources. 

2.2.2. Limiting deptlis of tine catalogs and flux uncertainties 

The noise in the Herschel catalogs results from the combined ef- 
fects of (1) instrumental effects + photon noise, (2) background 
fluctuations due to the presence of sources below the detection 
threshold (photometric confusion noise, see Dole et al. 2004), (3) 
blending due to neighboring sources, above the detection thresh- 
old (source density contribution to the confusion noise). In both 
PACS and SPIRE images (except at lOOyum in GOODS-N), the 



depths of the GOODS-Herschel observations are always lim- 
ited by confusion, i.e., (2) and (3) are always stronger than (1). 
Global confusion limits have been determined for PACS (Berta 
et al. 2011) and SPIRE (Nguyen et al. 2010). However these 
global definitions assume no a priori knowledge on the local 
projected densities of sources, as if e.g., 500 //m sources were 
distributed in an independent manner with respect to shorter 
wavelengths such as the 250 and 350yum ones, or even down 
to 24 jum. Moreover, the flux limit associated to source blending, 
(3), is often artificially set to be the flux density above which 
10% of the sources are blended, even though statistical studies, 
such as the present one, could afford higher fractions as long as 
the photometric uncertainty is well controlled. Actual observa- 
tions instead demonstrate that shorter wavelengths do provide 
a good proxy for the density field of longer wavelengths (see 
Fig.[T]i. Hence we define the 3cr (or 5cr) sensitivity limits of the 
GOODS-Herschel catalogs as the flux densities above which at 
least 68 % of the sources can be extracted with a photometric 
accuracy better than 33 % on the basis of Monte Carlo simula- 
tions and we use the positions of 24 pm sources as priors to ex- 
tract sources from PSF fitting. Individual sources are attributed a 
"clean" flag depending on the underlying density field as defined 
in Sect. 123] 

Flux uncertainties were derived in two independent ways. 
First we added artificial sources into the real Herschel im- 
ages and applied the source extraction procedure. This process 
was repeated a large number of times (Monte Carlo - MC - sim- 
ulations). Second we measured the local noise level at the 
position of each source on the residual images produced after 
subtracting sources detected above the detection threshold. The 
first technique gives a noise level for a given flux density aver- 
aged over the whole map, while the second one provides a local 
noise estimate. In the MC simulations, we define the 3cr (or 5cr) 
sensitivity limits in all bands as the flux densities above which a 
photometric accuracy better than 33 % (or 20 %) is achieved for 
at least 68 % of the sources in the faintest flux density bin (as in 
Magnein et al. 2009, 2011). 

Technique ( /) provides a statistical noise level attributed for a 
given flux density which accounts for all three noise components 
but is independent of local variations of the noise. The histogram 
of the output - input flux densities of the MC simulations follows 
a Gaussian shape whose rms was used to define the typical lim- 
iting depths of the Herschel catalogs listed in Col. (3) of Table[T] 
AU GOODS-//erac/ze/ images (except the PACS-lOOyum image 
in GOODS-N) reach the 3cr confusion level, i.e., the flux den- 
sity for which the photometric accuracy is better than 33 % for 
at least 68 % of the sources is more than three times higher than 
the instrumental noise level. 

In technique only the noise components (1) and (2) are 
taken into account, since the objects participating in the third 
component (source blending) have been subtracted to produce 
the residual images. However, imperfect subtraction of sources, 
due to local blending, may inflate the local residuals in the maps 
after source subtraction. In the PACS images and catalogs, both 
techniques result in very similar noise levels. A statistical limit- 
ing depth was computed by convolving the residual images with 
the PACS beam at each wavelength and measuring the rms of 
the distribution of individual pixels. This method resulted in the 
same depths as in technique (i) and listed in Col. (3) of Table [T] 
Instead, for the SPIRE data, local noise estimates in the residual 
maps were found to be systematically lower than those measured 
with technique On average, sources with a SPIRE flux den- 
sity corresponding to the detection threshold of 3 cr in the MC 
simulations are found to present a local signal-to-noise ratio of 
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Notes. Column definitions: Col.(l) central wavelength of the passband; Col. (2) full width half maximum (FWHM) of the point spread function 
(PSF) in the passband. In the shortest bands, the FWHM is limited by the under-sampling of the PSF; Col. (3) depth of the image at that wavelength, 
i.e., flux density of the faintest sources of the catalog. The depths listed in the table correspond to the 3tr limit. Due to local noise variations in the 
maps, some sources with slightly fainter flux densities may lie above this signal-to-noise threshold; Col. (4) total number of point sources above the 
3-(T limit. For the Herschel PACS and SPIRE passbands, we list the number of sources identified using a PSF-fitting based on Spitzer-MWS 24yum 
prior positions, themselves resulting from Spitzer-IKAC 3.6yum priors; Col. (5) number of sources identified with an optical counterpart having a 
spectroscopic redshift; Col. (6) fraction of sources with an optical counterpart having a spectroscopic redshift; Col. (7) fraction of sources with an 
optical counterpart having either a spectroscopic or a photometric redshift; Col. (8) number of sources used in the present study, i.e., sources which 
are both "clean" (not polluted by bright neighbors as discussed in Sect. I2.3t and for which a redshift either photometric or spectroscopic was 
measured; Col. (9) fraction of the sources listed in column (8) for which a spectroscopic redshift was determined. Cols. (10) to (16) for GOODS- 
south are defined as Cols.(3) to (9) for GOODS-north. 

The characteristics of the GOODS-Herschel PACS 100, 160 yum and SPIRE 250, 350 and 500 yum images and catalogs are given in the bottom 
part of the table. The SPIRE images extend over 30'x30' but we restricted the present analysis, hence also the number of sources given in the 
table, to the same field as the one covered with the PACS bands and to the sources with a 24 yum counterpart. The upper part of the table lists the 
characteristics of the Spitzer IRAC 3.6, 4.5, 5.8, 8yum, IRS peakup array 16;um and MIPS 24, 70yum images and catalogs. 

The depth of the SPIRE catalogs given here applies only to the sub-sample of "clean" galaxies, located in isolated areas as probed by the density 
maps obtained from the shorter wavelengths starting at 24 pm (see Sect. 12. 3b . This explains why they are much lower than the statistical confusion 
limits as listed by Nguyen et al. (2010) of 29, 31 and 34 mJy (5cr^(,„f confusion limits). 



5 in the residual maps. For SPIRE sources, this implies that we 
consider only sources above the 5cr limit in the residual maps, 
to be consistent with the 3cr limit resulting from the MC simula- 
tions. 

Due to local noise variations in the maps, there can be small 
numbers of sources with flux densities slightly fainter than the 
nominal detection limits, which explains the presence of sources 
below the horizontal lines in Fig.|5] 

2.3. Local confusion limit and "clean index" 

The main source of uncertainty, in the SPIRE images in particu- 
lar, comes from the high source density relative to the beam size, 
i.e., the so-called confusion limit (see Condon 1974). Assuming 
that this limit applies equally at all positions of the sky, Nguyen 
et al. (2010) estimated that the floor below which SPIRE sources 
may not be extracted is ~ 30 mJy, corresponding to 5<Tconf con- 
fusion limits of 29, 31 and 34 mJy/beam for beams of 18.1", 
24.9"and 36.6" FWHM at 250, 350 and 500 yum respectively. 

However, this "global confusion limit" is defined assuming 
no a priori knowledge on the projected density map of the un- 
derlying galaxy population. If one instead assumes that shorter 
wavelengths, at a higher spatial resolution, can be used to define 
the local galaxy density at a given galaxy position, then a "lo- 
cal confusion limit" can be defined. In practice, this means that 
not all SPIRE sources are located at a place where several bright 



PACS or MIPS-24yum fall in the SPIRE beam. Following this 
recipe, Hwang et al. (2010a) defined a "clean index" that was 
attributed to all individual Herschel detections under the fol- 
lowing conditions: a 500 //m source is flagged as "clean" if its 
24 jjm prior has at most one bright neighbor in the Spitzer-MlPS 
24/im band (where "bright" means an F24>50% of the central 
24;um source) within 20" (l.lxFWHM of Herschel at 250 yum) 
and no bright neighbor in each one of the shorter Herschel pass- 
bands, i.e., at 100, 160, 250, 350 and 500 yum within l.lxFWHM 
of Herschel in these passbands (see Table [TJ. As a result, we 
only kept 1 1 clean sources at 500 yum for which we consider that 
the photometry is reliable. The criterion becomes less critical 
for the shorter bands, since we only consider the presence of 
bright neighbors at shorter wavelengths. As a result, the num- 
ber of 350//m detections is an order of magnitude larger than 
at 500 fim. This "local confusion limit" was empirically defined 
after visually inspecting the data for all individual sources but a 
more detailed investigation of this quality flag using simulations 
of the actual GOODS sources both spatially and in redshift con- 
firms its robustness (Leiton et al. 2011, in prep.). For galaxies 
for which this "clean index" condition is not met in some bands, 
unphysical jumps in the IR SED are observed. This may lead to 
wrong estimates of the dust temperature for example, systemat- 
ically shifting it to colder values, since source blending aff'ects 
preferentially the longest wavelengths. 



6 



D. Elbaz et al.: GOODS-Herschel: an infrared main sequence for star-forming galaxies 




redshift redshift redshift 






Fig. 5. Distribution of the 'clean' GOODS-Herschel and Spitzer-MlPS 24 fim flux densities as a function of redshift (spectroscopic- 
filled and photometric-open symbols). Non "clean" galaxies (see Sect. I2.3l l are not represented. The average 3-cr depths of the 
surveys listed in Table [T] are shown with horizontal lines (the dashed lines at 100 and 160;um represent the shallower GOODS-N 
depths). Due to local noise variations, some 3-cr sources can be detected below these typical limiting flux densities. Sources above 
z~3 are not represented due to redshift uncertainties. 



With the sensitivity limits of GOODS-//ersc/ze/, Fig. |4] 
shows that below a redshift of z~3 the shortest wavelengths are 
always deeper than the longest ones, hence one can take advan- 
tage of these higher resolution images to better constrain the 
confusion limit at local, instead of global, scales. Moreover, we 
note that the fluxes in e.g., the SPIRE bands are not independent 
of those measured in the 24 //m and 100 /im passbands. They 
even follow a tight correlation (see Elbaz et al. 2010 and the 
present analysis), again in the redshift range of interest here, i.e., 
z~0-2.5. Hence it is possible to map the density of IR sources 
and to flag sources in relatively isolated areas, with respect to 
similarly bright or brighter IR sources. If the "clean index" did 
not reject efficiently problematic measurements, this would re- 
sult in an increase of the dispersion in the figures presented in 
this paper Since we will show that these dispersions are quite 
small already, if this effect was corrected, it would only rein- 
force our results. Typically, half of the Herschel sources detected 
at /l>160/im survive this criterion (see Table[T). 

3. High- and low-redshift galaxy samples 

3.1. The GOODS-Herschel ga\axy sample 

Both GOODS fields have been subject to intensive follow-up 
campaigns, resulting in a spectroscopic redshift completeness 
greater than 70 % for the Herschel sources (Table [TJ. We use 
a compilation of 3630 and 3018 spectroscopic redshifts for 
GOODS-N (Cohen et al. 2000, Wirth et al. 2004, Barger, Cowie 
& Wang 2008, and Stem et al. in prep.) and GOODS-S (Le 
Fevre et al. 2004, MignoU et al. 2005, Vanzella et al. 2008, 
Popesso et al. 2009, Balestra et al. 2010, Silverman et al. 2010, 
and Xia et al. 2010) respectively. Photometric redshifts and stel- 
lar masses are computed in both fields from U-band to IRAC 
4.5 /im photometric data using Z-PEG (Le Borgne & Rocca- 
Volmerange 2002). The templates used for both photometric red- 



shifts and stellar mass estimates are determined from PEGASE.2 
(Fioc & Rocca-Volmerange 1999) are were produced using nine 
scenarios for the star formation history (see Le Borgne & Rocca- 
Volmerange 2002) with various star-formation efficiencies and 
infall timescales, ranging from a pure starburst to an almost con- 
tinuous star-formation rate, aged between 1 Myr and 13 Gyr (200 
ages). There is no constraint on the formation redshift. The tem- 
plates are required to be younger than the age of the Universe at 
any redshift. 

The redshift distributions of the sources individually de- 
tected in each of the Herschel bands, as well as at 24 yum with 
Spitzer, are presented in Fig. |5]for both fields. This illustrates the 
relative power of these bands to detect sources as a function of 
redshift. While the 500 fim band samples sources at all redshifts 
from z=0 to 4, it only provides a handful of objects: 24 galaxies 
in total within the 10'xl5' size of the GOODS-N field, with only 
1 1 flagged as clean, 73 % of which have a spectroscopic redshift 
determination. In comparison, more than a thousand sources are 
detected in the 100//m band, the vast majority being flagged as 
clean and 72 % having a spectroscopic redshift. 

Table [T] also lists the characteristics of the other IR catalogs 
that we use in the present study. The GOODS Spitzer IRAC cata- 
logs were created using SExtractor (Bertin & Arnouts 1996), de- 
tecting sources in a weighted combination of the 3.6 and 4.5 yum 
images, with matched-aperture photometry in the four IRAC 
bands, using appropriate aperture corrections to total flux. The 
Spitzer 24yum and lOjjm catalogs (Magnelli et al. 2011) use 
data from the Spitzer GOODS and FIDEL programs (PI: M. 
Dickinson). Sources detected in the IRAC images are used as 
priors to extract the 24 fim fluxes, and then in turn a subset of 
those 24 jjm sources are used as priors to extract fluxes at 70 yum. 
The 16 fim data comes from Spitzer IRS peak-up array imaging 
(Teplitz et al. 2011); here again, 16/im catalog fluxes are ex- 
tracted using IRAC priors. In this study, we make particular use 
of the Spitzer data to quantify the redshift dependence of the IR 
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SEDs while minimizing mid-infrared k-corrections by measur- 
ing the rest-frame 8 fim emission of galaxies at z~0, 1 and 2 from 
their observed fluxes in the IRAC-8 jjm, IRS-16yum and MIPS- 
24//m passbands. Table[T]also gives the spectroscopic {%Zspec), 
and photometric + spectroscopic i%Zaii) completeness of the IR 
catalogs from 3.6 to 500 //m within the fiducial GOODS area. As 
noted previously, the SPIRE images of GOODS-N cover a wider 
field, but here we do not count the sources detected outside the 
regular GOODS area. 

Known AGN were excluded from the sample and will be dis- 
cussed separately in Section |8] X-ray/optical AGN were iden- 
tified from one of the following criteria: Lx[0. 5-8.0 keV] > 
3x10"^^ ergs s"', a hardness ratio (ratio of the counts in the 2-8 
keV to 0.5-2 keV passbands) higher than 0.8, Nh >10^^ cm"^, or 
broad/high-ionization AGN emission lines (Bauer et al. 2004). 
We also excluded power-law AGN, i.e., galaxies showing a ris- 
ing continuum emission in the IRAC bands due to hot dust radi- 
ation (see definition in Sect.[8]l. 



3.2. Total infrared luminosities 

Total IR luminosities, L»"'*'=*"=', for GOODS-Herschel galaxies 
were determined by allowing the normalization of the CEOl tem- 
plate SEDs to vary and choosing the one that minimizes the;^'^ fit 
to the Herschel measured flux densities. At the highest redshifts 
considered in the present analysis (z k 2.5), the Herschel 100 yum 
passband samples rest-frame mid-IR wavelengths. Hence, to 
avoid mixing galaxies with and without direct far-IR detections, 
we require at least one photometric measurement at wavelengths 
longer than 30;um in the rest-frame. This excludes a few high 
redshift galaxies detected only at lOOyum. Total IR luminosities, 
^iteischei^ were integrated from 8 to 1000/zm on the best-fitting 
normalized CEO I SED. When only one or two Herschel mea- 
surements are available above 30 yum, the degeneracy of the fit 
being large, we use the standard CEOl technique, i.e., we use 
the SED with the closest luminosity from the CEOl library with- 
out allowing any renormalization. 

In order to quantify the impact of the choice of a given set of 
SEDs to fit the Herschel measurements and determine Lj^'"^*'^'"^', 
we have repeated the same exercise with another SED library 
from Dale & Helou (2002, DH02). The ratio of the LH'^rschei ^^j. 
ues derived with one or the other family of SEDs has a median of 
1 and a dispersion of 12 %-rms. The uncertainty in the determi- 
nation of Lj^"schei jg therefore dominated by the actual error bars 
on the Herschel flux measurements rather than by the choice of 
the SED library. In order to account for the latter source of un- 
certainty, we have generated a series of 100 realizations of the 
Herschel flux measurements assuming a Gaussian distribution 
within their error bars and determined 100 values of /.j^'^rschei 

asso- 



Table 2. Number of galaxies and total IR luminosity range of the 
local galaxy samples. 



IR 



fitting those realizations independently. The final L\ 
ciated to a given galaxy is the median of the 100 Monte Carlo 
estimates and its error bar is the rms around the median. This 
procedure was repeated for each individual galaxy. 

Since we will compare the distant GOODS-Herschel galax- 
ies to a reference sample of local galaxies for which Lj°' is es- 
timated from IRAS measurements alone, as a consistency check 
we computed the total IR luminosity that we would obtain for 
the GOODS-Herschel galaxies if we had used Eq.[T](taken from 
Sanders & Mirabel 1996), 



Local 
samples 


(2) 
Nb 


(3) 

Rradio 


(4) 
FEE 


<10 


logio 
10-11 


(iiR/Lo) 
11 - 12 


>12 


ISO 


150 


11 





36 


56 


14 


45 


AKARI 


287 


47 





63 


164 


55 


9 


Spitzer 


211 


58 


211 





44 


154 


13 


Total 


648 


116 


211 


99 


264 


223 


67 



where FIR 



turn 



m] [l.8 X lO-^\FIR[Wm-^])\ /3.826 x 10' 

l3.4SFi2fjm + 5.16F25//m + '2.5SF(,()^,„ + Fnt)^,„ 



,26 



Notes. Column (2) total number of objects for each local sample; 
Col. (3) number of galaxies with a radio continuum (1.4 GHz) size esti- 
mate; Col. (4) number of galaxies for which a fraction of extended emis- 
sion (FEE) was measured (Diaz-Santos et al. 2010), i.e., fraction of the 
mid-IR continuum at 13.2;um more extended than the Spitzer/IKS res- 
olution of 3.6 " (see Sect.|5}. 



as a proxy for the derivation of the 8 - 1000 fim luminosity, in- 
stead of the actual integral over the IR SED . The IRAS flux densi- 
ties Fn^lm, F2Siim, Fa)^,„ and Fm,,m in Eq.fflare in Jy. Both tech- 
niques give equivalent total IR luminosities within 5 %, hence 
again the dominant cause of discrepancy in the comparison is 
related to flux uncertainties. 



3.3. Local galaxy reference sample 

The local galaxy reference sample that we use in this paper con- 
sists of galaxies detected with the Infrared Space Observatory 
(ISO), AKARI, and Spitzer. Their rest-frame 8 pm luminosities 
and total IR luminosities are compared to those of the GOODS- 
Herschel galaxies. Galaxies with direct IRAC-8 pm measure- 
ments from Spitzer are supplemented with galaxies with ISO 
6.75yum and AKARI 9pm photometry, for which pseudo-IRAC 
8 pm luminosities, Lg, were computed using the IR SED of M82 
(Forster Schreiber et al. 2001, Elbaz et al. 2002). The ISO and 
AKARI samples span a wide range of relatively low luminosity 
galaxies, together with a sample of ULIRGs, while the Spitzer 
sample contains a quite complete sample of local LIRGs (see 
TableEll. 

3.3.1 . Local ISO galaxy sample 

The mid-IR luminosities of this sample of 150 galaxies de- 
scribed in CEOl and Elbaz et al. (2002) were obtained from 
measurements taken with ISO. The sample includes 1 10 galaxies 
closer than 300 Mpc and spanning a wide range of mid-IR lumi- 
nosities estimated from ISOCAM-LW2 (5-8.5//m, centered at 
6.75 pm) and 41 ULIRGs, at distances 80 to 900 Mpc, with mid- 
IR luminosities determined with the PHOT-S spectrograph of 
ISOPHOT (Rigopoulou et al. 1999). We refer to CEOl for a dis- 
cussion of the conversion of the PHOT-S spectra into broadband 
luminosities equivalent to the LW2 filter Pseudo-IRAC 8 pm lu- 
minosities, Ls, were estimated by first convolving the ISOCAM 
CVF spectrum of M82 (Forster Schreiber et al. 2001, Elbaz et 
al. 2002) to the ISOCAM-LW2 and IRAC-8yum bandpasses and 
then normalizing the resulting luminosities to the observed lu- 
minosity for each of the 150 galaxies, in order to derive their Lg. 
Since both filters are wide and largely overlapping, the conver- 
sion depends very little on the exact shape of the spectrum used 
for the conversion and we checked that indeed using the CEOl 
SEDs (for example) instead of that for M82 would make negli- 
gible differences with respect to the actual dispersion of galaxies 
antheL|°'- 



.^y w ^jjj Lg diagram. Total IR luminosities, L|°', were derived 

from the four IRAS band measurements using Eq. [T] 
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3.3.2. Local AKARI galaxy sample 

Galaxies with mid-infrared (9 fim) measurements from AKARI 
were cross-matched with the IRAS Faint Sources Catalog ver 
2 (FSC-2; Moshir, Kopman & Conrow 1992) and with spectro- 
scopic redshifts from the Sloan Digital Sky Survey Data Release 

7 (SDSS DR7; Abazajian et al. 2009) supplemented by a pho- 
tometric sample of galaxies with redshifts available in the lit- 
erature (Hwang et al. 2010b). For both IRAS and AKARI, we 
consider only the sources with reliable flux densitiefl A total of 
287 galaxies have 9 fim flux densities from the A/iTA/JZ/Infrared 
Camera (IRC, Onaka et al. 2007) Point Source Catalog (PSC ver. 
1.0, Ishihara et al. 2010) reaching a detection limit of 50 mJy 
(5cr) with a uniform distribution over the whole sky and closer 
than -450 Mpc (z<0.1). As in Sect. lTlTl pseudo-IRAC 8 yum lu- 
minosities, Lg, were computed by convolving the ISOCAM CVF 
spectrum of M82 with the AKARI-IRC 9 yum bandpass to esti- 
mate the conversion factor between the IRC-9 //m and IRAC- 

8 jum luminosities assuming the same IR SED for all galaxies. 
The effective wavelength of the AKARI 9 jum passband is 8.6 yum 
(Ishihara et al. 2010), not far from that of the IRAC-8 fim filter 
(7.9jum, Fazio et al. 2004). 

Total IR luminosities were computed from the four IRAS 
bands using Eq. [1] The IRC-9 yum measurements were not used 
in the computation of Lj°'. Far-IR measurements were supple- 
mented with the AKARIfFca'-lnfmred Surveyor (FIS; Kawada et 
al. 2007) all-sky sui-vey Bright Source Catalogue (BSC ver 1 .C0) 
that contains 427 071 sources, with measured flux densities at 
65, 90, 140 and 160 jim. We used the supplementary far-IR mea- 
surements for 16 % of the sample for which there is no 12 /im 
nor 25 fim reliable measurement from IRAS. We checked the 
consistency of these IR estimates from AKARI with those ob- 
tained from IRAS alone and found that AKARI luminosities were 
systematically lower by 10 %. We corrected those 16 % galaxies 
by this factor. 

3.3.3. Local Spitzer galaxy sample 

A sample of 202 IRAS sources, consisting of 291 individ- 
ual galaxies (some blended at IRAS resolution), were observed 
with the IR spectrograph (IRS) on-board Spitzer as part of the 
Great Observatories All-sky LIRG Survey project (GOALS; 
Armus et al. 2009). The sources were drawn from the IRAS 
Revised Bright Galaxy Sample (RBGS; Sanders et al. 2003) 
and represent a complete sub-sample of systems iz < 0.088) 
with IR luminosities originally defined to be in the range 
of lO"L0<LiR<lO'^L0. The GOALS sample includes 200 
LIRGs and 22 ULIRGs. The total IR luminosities of the sys- 
tems were derived using their IRAS measurements and Eq. 1 (see 
Armus et al. 2009 for further details on this calculation). 

Using the spectral images obtained with the short-low mod- 
ule of IRS, Dfaz-Santos et al. (2010) measured the spatial extent 
of the light radiated in the mid-IR continuum at 13.2/vm of a 
sub-sample of 21 1 individual galaxies (closer than 350 Mpc) for 
which data were available at the time of publication and sources 
could be detected. We use these size estimates in our analysis re- 
garding the link between star formation compactness and the //?8 
ratio. This fraction of extended emission (FEE) is directly re- 
lated to the spatial distribution of the star formation regions and 
presents the advantage of being measured in a wavelength range 

' Flux quality flags are either "high" or "moderate" for IRAS sources 
and "high" for AKARI sources 

2 http://www.ir.isas.jaxa.jp/AKARI/Observation/PSC/Public/'RN/AK. 
FIS_BSC_Vl.RN.pdf 



not affected by the presence/absence of emission lines such as 
R\Hs. For the multiple systems unresolved by IRAS, Dfaz- 
Santos et al. (2010) distributed the total IR luminosity between 
galaxies proportionally to their SpitzerjMS'S-lA pva fluxes. Due 
to this redistribution of the luminosity, there are now 44 galaxies 
with IR luminosities less than 10" Lq in our sample. Added to 
these normal star-forming galaxies, the present sample finally in- 
cludes 154 LIRGs and 13 ULIRGs (with 10'2<L[°*/Lo<4xl0'2). 
IRAC-8 jum luminosities for these galaxies are from Mazzarella 
et al. (in prep). Stellar masses were derived by cross-matching 
the GOALS sample with 2MASS and converting the Ks lumi- 
nosities into stellar masses (excluding remnants) using a us- 
ing a mass-to-light ratio M^ILks-Q-I Mq/L^:,© computed from 
PEGASE 2 (Fioc & Rocca-Volmerange 1997, 1999) assuming a 
Salpeter IMF and an age of 12 Gyr 

4. Universality of 7/^8 (=LiR/Lg): an IR main sequence 

4.1. The mid-infrared excess problem 

Before the launch of Herschel, the derivation of Lj°', hence 
also of the SFR, of distant galaxies had to rely on extrapola- 
tions from either mid-IR or sub-mm photometry. While there 
are many reasons why extrapolations from the mid-IR could be 
wrong (evolution in metallicity, geometry of star formation re- 
gions, evolution of the relative contributions of broad emission 
lines and continuum), it was instead found that they work rela- 
tively well up to z~1.5. Using shallower Herschel data than the 
present study, Elbaz et al. (2010) compared L™, estimated from 
Herschel PACS and SPIRE, to Lj^ - the total IR luminosity ex- 
trapolated from the observed Spitzer mid-IR 24 jim flux density - 
and found that they agreed within a dispersion of only 0.15 dex. 
The CEOl technique used to extrapolate L^^ attributes a single 
IR SED per total IR luminosity. Hence a given 24/im flux den- 
sity is attributed the L|°' of the SED that would yield the same 
flux 24 yum flux density at that redshift. 

Stacking Spitzer MIPS-70yum measurements at prior posi- 
tions defined by 24yum sources in specific redshift intervals, 
MagneUi et al. (2009) found that the rest-frame 24 yum/( 1 -Hz) 
and 70yum/(l+z) luminosities were perfectly consistent with 
those derived using the CEOl technique for galaxies at z<1.3. 
Although the 70 yum passband probes the mid-IR regime for red- 
shifts z>0.8, it presents the advantage of sampling the contin- 
uum IR emission of distant galaxies without being affected by 
the potentially uncertain contribution of PAHs, contrary to that 
at 24 yum. At z> 1 .5 however, extrapolations from 24 fim measure- 
ments using local SED templates were found to systematically 
overestimate the 70 /im measurements (Magnelli et al. 2011). 
This mid-IR excess, first identified by comparing Lj^ with radio, 
MIPS-70//m and 160 yum stacking (Daddi et al. 2007a, Papovich 
et al. 2007, Magnelli et al. 2011) has recently been confirmed 
with Herschel by Nordon et al. (2010) on a small sample of z~2 
galaxies detected with PACS and by stacking PACS images on 
24fj.m priors (Elbaz et al. 2010, Nordon et al. 2010). 

Here, thanks to the unique depth of the GOODS-Herschel 
images, we are able to compare Lj^ to for a much larger num- 
ber of galaxies than in Elbaz et al. (2010) and, more importantly, 
for direct detections at z>1.5. In the left-hand part of Fig. |6] we 
show that the mid-IR excess problem is not artificially produced 
by imperfections that could result from the indirect stacking 
measurements, but instead takes place for individually detected 
galaxies at z> 1 .5 and at high 24 pm flux densities, coiTespond- 
iig to L^^>10'^ Lq. Although known AGN were not included 
in the sample, unknown AGN may still remain. Indeed it has 
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Fig. 6. Left: Comparison of L}^ (8-1000 yum) as directly measured from Herschel (LJ^"™"') with the value extrapolated from lAjim 
(L^^'*^^"') using the CEOl technique. Only "clean" galaxies are represented (as defined in Sect. 13.1b . Galaxies with spectroscopic 
and photometric redshifts (from both GOODS-north and south) are marked with filled and open symbols respectively. Colors range 
from black (z~0) to orange (z~2.5), passing through green (z~l) and red (z~2). The wavelength range sampled by the MIPS- 
24/im passband is shown in orange at the top of the figure where it is compared to the redshifted SED of M82. The dashed line 
in the left-hand side panels is the one-to-one correlation. The sliding median and 16th and 84th percentiles of the distribution are 
shown with white dots connected with a solid line and grey zone respectively. The bottom panel shows the ratio of the actual 
over extrapolated total IR luminosity. Right: Comparison of /.j^'"'*^' with (rest-frame 8 im\ broadband) for "clean" galaxies. 
The observed bandpasses used to estimate Lg are illustrated in the top of the figure and compared to the redshifted SED of M82. 
The sliding scale of the median and 68 % dispersion around it is shown with a grey zone which is fitted by the solid and dashed 
lines: //?8=4.9 [-2.2,-1-2.9]. Stacked measurements combined with detections weighted by number of objects per luminosity bin are 
represented by large yellow open triangles (GOODS-south: upside down, GOODS-north: upward). The bottom panel shows the 
IR% (=Lj^"^'"'Vi'8) ratio which is found to remain constant with luminosity and redshift. 



been proposed that the mid-IR excess problem could be due to 
the presence of unidentified AGN affected by strong extinction, 
possibly Compton thick (Daddi et al. 2007b, see also Papovich 
et al. 2007). At these high redshifts, the re-processed radiation of 
a buried AGN may dominate the mid-IR light measured in the 
24 yum passband, while the far-IR emission probed by Herschel 
would be dominated by dust-reprocessed stellar light. Indeed, 
studies of local dusty AGN have demonstrated that their contri- 
bution to the IR emission of a galaxy drops rapidly above 20 /im 
in the rest-frame (Netzer et al. 2007). However, this explanation 
for the mid-IR excess problem was recently called into ques- 
tion by mid-IR spectroscopy of z~2 galaxies obtained using the 
Spitzer IRS spectrograph showing the presence of strong PAH 
emission lines where one would expect hot dust continuum emis- 
sion to dominate if this regime were dominated by a buried AGN 
(Murphy et al. 2009, Fadda et al. 2010) and by deeper Chandra 
observations (Alexander et al. 201 1). 



4.2. Resolving the mid-IR excess problem: universality ofIR8 

We have seen that extrapolations of Lj°' from 24 /im measure- 
ments using the CEOl technique fail at z>1.5. We also find that 
using the same technique with another set of template SEDs, 
such as the DH02 ones, fails in a similar way. 

We wish to test the main hypothesis on which the CEOl tech- 
nique relies, namely, that IR SEDs do not evolve with redshift. If 
that was the case, then a single SED could be used to derive the 
Lj°' of any galaxy whatever the rest-frame wavelength probed, 
as long as it falls in the dust reprocessed stellar light wavelength 
range. Indeed, local galaxies are observed to follow tight corre- 
lations between their mid-IR luminosities at 6.75, 12, 15, 25 //m 
and L',°' (see CEOl, Elbaz et al. 2002) as well as with their SFR 
as derived from the Paff line (Calzetti et al. 2007) for the Spitzer 
passbands at 8 and 24 /vm. This technique fails at z>1.5, which 
has until now been interpreted as evidence that distant IR SEDs 
are different from local ones. However, in order to properly test 
the redshift evolution of the IR SEDs, it is necessary to compare 
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measurements in the same wavelength range for galaxies at all 
redshifts. For that purpose, we now compute the same rest-frame 
mid-lR luminosity, (-vLy[^ /um]), defined as the luminos- 
ity that would be measured in the lRAC-8 fim passband in the 
rest-frame. We choose this particular wavelength range because 
it can be computed from z~0 to 2.5 with minimum extrapola- 
tions using the lRAC-8 fim filter for nearby galaxies (z<0.5), the 
IRS-16/im peak-up array for intermediate redshifts around z~l 
(0.5<z<1.5) and the MIPS-24yum passband at z~2 (1.5<z<2.5). 
Even in these conditions, small k-corrections need to be applied 
in order to calculate for the same rest-frame passband. This 
was done using the mid-IR SED of M82 for all galaxies. We 
verified that using other SEDs, such as the CEOl or DH02 tem- 
plates, would alter Lg by factors that are small when compared 
with the dispersion of the observed - Lg relation. The re- 
sults are shown in the right-hand part of Fig. |6] Surprisingly, 
when plotting galaxies at all redshifts and luminosities in the 
same wavelength range, we no longer see a discrepancy between 
galaxies above and below z~1.5. The sliding median of the IRS 
ratio, defined as IRS=L}^/Ls, - illustrated by white points con- 
nected with a solid grey line in the right-hand part of Fig. |6]- 
remains flat and equal to IRdi-4.9 [-2.2,-1-2.9] (solid and dashed 
lines in Fig.|6]-right) from L8=10^ to 5x10" Lq or equivalently 
from Lj°'=5xlO'' to 3x10'^ Lq. The 68 % dispersion around the 
median is only +0.2 dex. 

In order to test possible selection effects on the galaxies used 
to determine the //?8 ratio, we combined Herschel detections 
with stacked measurements on lAjiva prior positions. This was 
done by defining intervals of luminosity in Lg, e-g-, from the 
16 jum band for sources around z~ 1 or 24 jum for sources around 
z~2. In a given Lg interval, we determined the median of the Lj°' 
obtained for detections on one hand (white dots connected with a 
solid line in the right-hand part of Fig. |6) and on the other hand 
measured average PACS 100 yum and 160 yum flux densities for 
the sources with no Herschel detection by stacking sub-images 
of 60" on a side at their 24 jim prior positions. These sub-images 
were extracted from the residual images to avoid contamination 
by detections. The average stacked PACS-lOOjum and 160 /im 
flux densities were converted into total IR luminosities using the 
CEOl library of template SEDs, selected based on luminosity at 
the median redshift of the galaxies in that Lg luminosity interval. 
We found no systematic difference when deriving Lj^' from the 
PACS 100 jum or 160 jum data when using the CEOl templates 
for the extrapolation (see also Elbaz et al. 2010) . Both PACS 
bands gave consistent values for Lj°'. The two values obtained 
for Lj°' from detected and stacked undetected sources were then 
combined according to a weight depending on the number of 
sources in each group within this Lg interval and on the signal- 
to-noise ratio of these measurements (quadratically), in order to 
avoid giving the same weight to both measurements if they have 
the same number of sources but very different S/N ratios. The 
resulting - Lg relation is shown with yellow open triangles 
separately for each GOODS field. Since the 100 /im and 160 /im 
gave similar results, we only present in the right-hand part of 
Fig.|6]the result obtained from the 100 /im band. Again, the typ- 
ical IR% ratio appears to be flat, independent of both luminos- 
ity and redshift. The range of luminosities probed by GOODS- 
Herschel varies as a function of redshift as shown in the upper 
panel of Fig. |7] where we represent the distribution of total IR 
luminosities measured with Herschel as a function of redshift 
for the galaxies classified as "clean" (Sect.|2]i. This is due to the 
combination of limited volume at low redshifts - limiting the 
ability to detect rare luminous objects - and depth at high red- 
shifts - limiting the ability to detect distant low luminosity ob- 
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Fig. 7. Upper panel: Distribution of total IR luminosities 
(Lj^™'""') of the GOODS-Herschel galaxies classified as "clean" 
as a function of redshift. The solid and dashed red lines are the 
detection limits of the GOODS-S and GOODS-N images re- 
spectively. Spectroscopic and photometric redshifts are shown 
with filled and open dots respectively. Bottom panel: //?8 ratio 
(=LJ^"'^'"°'/L8) as a function of redshift. The solid and dashed 
horizontal black lines are the median and 16th and 84th per- 
centiles of the distribution (Eq. |3]l, i.e., IR8^4.9 [-2.2, H-2.9]. 
The solid and dashed red lines show the detection limits of the 
GOODS-S and GOODS-N images above z=1.5. The sliding 
median of the sources detected by Herschel is shown with black 
open circles connected with a solid line. A weighted combina- 
tion of detections with stacked measurements (as in Fig. |6] and 
as described in Sect. 14.2b is shown with open yellow triangles 
(both fields combined). 



jects. In the bottom panel of Fig. |7] we show the redshift evolu- 
tion of the IRE ratio. It is flat up to z~2 and then, due to the shal- 
lower detection limit of Herschel compared to Spitzer-24 //m, it 
is slightly larger than the typical value, since only galaxies with 
high L'^/Ls can be detected by Herschel. 

Hence, we do not see a mid-lR excess when comparing sys- 
tematically Lj°' to 8/im rest-frame data. In particular, if AGN 
were playing a more important role at z> 1 .5 than at lower red- 
shifts, we would expect to see a change in //?8 at this redshift cut- 
off contrary to what is actually observed. The cause for the mid- 
IR discrepancy is therefore not specific to galaxies at z > 1.5, but 
is instead due to the templates used to represent z ~ ULlRGs. 
Locally, galaxies with > 10'^ Lq are very rare, most prob- 
ably because galaxies today are relatively gas-poor compared 
to those at high redshift. Moreover, they have infrared SEDs 
that are not typical of star-forming galaxies in general, includ- 
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Fig. 8. L5^'[IRAS] versus Lg for local galaxies including only the 
ISO sample of galaxies used to build the CEOl library of tem- 
plate SEDs and converted from 6.75 fim to 8 yum using the SED 
of M82. The Hght blue line shows the position of the CEOl SED 
templates, built to follow two power laws in the Lj°' - relation. 



ing those of most distant ULIRGs. The majority of high-redshift 
galaxies, even ultraluminous ones, share the same IR properties 
as do local, normal, star-forming galaxies with lower total lumi- 
nosities. Galaxies with SEDs like those of local ULIRGs do exist 
at high redshift, but they do not dominate high redshift ULIRGs 
by number as they do in the present day. 

4.3. Origin of the "mid-IR excess" discrepancy 

Fig. [8] shows the original Lj°'-L8 data that were used to build 
the CEOl library of template SEDs. The solid line in the figure 
shows the relation traced by the SED templates. Originally, the 
mid-IR luminosity was computed from the IS0CAM-LW2 filter 
at 6.75 fim, L^j, which we convert here to Lg using the SED of 
M82. The conversion was validated by a sub-sample of galaxies 
for which we have measurements with both IS0CAM-LW2 and 
IRAC-8 jim. While the ti-end followed by the CEOl templates is 
consistent with the GOODS-//ersc/ie/ galaxies below Lg ~10'" 
L0, there is a break above this luminosity threshold that was re- 
quired to fit the local ULIRGs in this diagram. 

In Fig. |9l-left, we supplement the original local ISO sam- 
ple with the 287 AKARI galaxies introduced in Sect. 13.3.21 With 
this larger sample, we see galaxies extending the low luminosity 
ti-end beyond the threshold of Lg ~10'° Lq, with a flat IRS ratio. 
This trend is similar to the one found for the GOODS-Herschel 
galaxies (background larger orange symbols as in Fig. |6) and the 
extended local sample is well contained within the 16th and 84th 



percentiles around the median of the GOODS-Herschel sample 
(solid and dashed lines in Fig.|6]l. 

The median of both samples are very similar (see Eqs. l3]2) . 

^^giocai^4g [-1.7, +6.4] (2) 
^^gGOODS-H.r.«M ^ 4 9 [-2.2, -h2.9] (3) 

Note, however, the large upper limit of the 68 % dispersion in 
Eq. |2] which is mainly due to the elevated //?8 values of the 
local ULIRGs, as seen in the left-hand panel of Fig.|9] The me- 
dians of both samples are shifted to higher values because of the 
asymmetric tails of galaxies with large values of IRS, as shown 
in the right-hand part of Fig. |9] where we compare the IRS dis- 
tribution for the local ISO+AKARI galaxies (upper panel) with 
that of the GOODS-Herschel galaxies (lower panel). Both dis- 
tributions present the same properties: they can be fitted by a 
Gaussian and a tail of high-//?8 values. The central values and 
widths cr of the Gaussian distributions are very similar for both 
samples (Eqs. l4l5l l, 

//?8'°""'(center Gaussian) = 3.9 [o-=1.25] (4) 
^^gGOODS-H.™M(^gjjjgj. Gaussian) = 4.0 [cr = 1.6] , (5) 

again reinforcing the interpretation that the distant galaxies be- 
have very similarly to local galaxies. If the IR SED of galaxies 
were diff'erent at low and high redshift, then one would not ex- 
pect them to have the same distributions in IRS. 

Hence we do not find evidence for diff'erent IR SEDs in dis- 
tant galaxies. Instead, we find that local and distant galaxies are 
both distributed in two quite well-defined regimes: a Gaussian 
distribution containing nearly 80 % of the galaxies, which share 
a universal IRS ratio of ~4, and a sub-population of ~20 % of 
galaxies with larger //?8 values. The exact proportion of this sub- 
population is not absolutely determined from this analysis, since 
it depends on the flux limit used to define the local reference 
sample, while the distant sample mixes together galaxies span- 
ning a large range of redshifts and luminosities. Nevertheless, 
the objects in the high-//?8 tail remain a minority at both low 
and high redshift compared with those in the Gaussian distribu- 
tion. 

In the following, we call the dominant population "main se- 
quence" galaxies, since they follow a Universal trend in L'j°'-Lg 
valid at all redshifts and luminosities. We also justify this choice 
in the next sections by showing that this population also follows 
a main sequence in SFR - M*, while galaxies with an excess 
IRS ratio systematically exhibit an excess sSFR (=SFR/M»). In 
the local sample, ULIRGs are clearly members of the second 
population whereas z~2 ULIRGs mostly belong to the Gaussian 
distribution, hence are main sequence galaxies. It is therefore the 
weight of both populations that has changed with time and that 
is at the origin of the mid-IR excess problem. The CEOl SED 
library, illustrated by a blue line in Figs. |8]and|9] reaches val- 
ues of IRS that are more than five times larger than the typical 
value for main sequence galaxies. This leads to an overestimate 
of L'j°' when the SED templates for local ULIRGs are used to 
extrapolate from 24 fim photometry for main sequence galaxies 
at z~2. Note, however, that it is not necessary to call for a new 
physics for the IR SED of these galaxies that would justify, e.g., 
stronger PAH equivalent widths, since most of the distant LIRGs 
and ULIRGs belong to the same main sequence as local normal 
star-forming galaxies. It is well-known that local (U)LIRGs are 
experiencing a starburst phase, with compact star formation re- 
gions, triggered in most cases by major mergers (see e.g., Armus 
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Fig. 9. Le/i;L|°'[IRAS] versus for local ISO and ^XA/?/ galaxies (filled blue dots). The GOODS-Herschel galaxies are shown 
in the background with lighter orange symbols (same points as in Fig. |6ll together with their median (black solid line) and 68 % 
dispersion (black dashed lines). The light blue line shows the locus traced by the CEOl SED library. Right: Histogram of the //?8 
ratios for the local galaxy sample (blue, upper panel) and GOODS-Herschel sample (orange, bottom panel). The solid curves show 
Gaussians fit to the distributions. The vertical grey lines indicate the median (solid) and 68 % dispersion (dashed) for the full samples 
(as in Eqs.[3]|2|. 



et al. 1987, Sanders et al. 1988, Murphy et al. 1996, Veilleux, 
Kim & Sanders 2002 for ULlRGs and Ishida 2004 for LIRGs). 
This leads us to the investigation of the role of compactness pre- 
sented in the next section. Indeed, if local ULIRGs are known 
to form stars in compact regions and are found to be atypical 
in terms of //?8, then it would be logical to expect that distant 
ULIRGs instead are less compact, perhaps as a result of their 
higher gas fractions. Note also that galaxies with an excess IRE 
ratio are found at all luminosities and redshifts and are not only 
a chai^acteristic of ULIRGs. 

5. //?8 as a tracer of star formation compactness 
and "starburstiness" in local galaxies 

The size and compactness of the star formation regions in galax- 
ies is a key parameter that can affect the IR SED of galaxies. 
Chanial et al. (2007) showed that the dust temperature (Tdust) 
estimated from the IRAS 60 over 100 jum flux ratio, R(60/100), 
is very sensitive to the spatial scale over which most of the IR 
light is produced. It is known that there is a rough correlation of 
R(60/100), hence Tdust, with L[°* (Soifer et al. 1987): locally, the 
most luminous galaxies are warmer This relation has recently 
been established with AKARI and Herschel in the local and dis- 
tant Universe (Hwang et al. 2010a). Locally, where galaxies can 
be spatially resolved in the far-IR or radio, Chanial et al. (2007) 
showed that the dispersion in the Lir - Tdust relation was signif- 
icantly reduced by replacing Lir by the IR surface brightness, 
Sir. We extend this analysis to the relation between this star for- 



mation compactness indicator, Sir, and IRE, the far-IR over mid- 
IR luminosity ratio. In the present study, the term "compactness" 
is used to refer to the overall size of the starburst and not to the 
local dumpiness of the various star formation regions, which we 
cannot measure in most cases. 

An extension of the Chanial et al. analysis to the brighter IR 
luminosity range of (U)LIRGS has become possible thanks to 
the work of Diaz-Santos et al. (2010). They used SpitzerflRS 
data to derive the fraction of extended emission of the mid- 
IR continuum of the GOALS galaxy sample (Sect. 13.3.31) at 
13.2 fim. 

5.1. Determination of tiie projected star formation density 

5.1 .1 . Radio/Far-IR projected surface brightness 

Due to the limited angular resolution of far-IR data, we first es- 
timate the sizes of star formation regions from radio imaging by 
cross-matching the local galaxy sample with existing radio con- 
tinuum surveys and then convert them into far-IR sizes using a 
correlation determined from a small sample of galaxies resolved 
in both wavelength domains as in Chanial et al. (2007). 

The IRAS-60jum and VLA radio continuum (RC, 20 cm) 
azimuthally averaged surface brightness profiles of a sample of 
22 nearby spiral galaxies was fitted by a combination of expo- 
nential and Gaussian functions by Mayya & Rengarayan (1997). 
The angular resolutions of the 60 yum and 20 cm maps used in 
that study was about 1', so we deconvolved their synthesized 
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Fig. 10. IRAS-60 /vm versus VLA 20 cm radio continuum half- 
light radius correlation. 

profiles by a 1' beam and derived the intrinsic half-light radii riR 
(at 60yL(m) and rRc- The half-light radii estimates at both wave- 
lengths are strongly correlated (Fig. [10); a logarithmic bisector 
fit to the data is given in Eq. |6l 

riR = (0.86 ± 0.05) trc (6) 

Hence in the following, we estimate the far-IR sizes of the star 
formation regions of our local galaxy sample from their radio 
continuum half-light radius using Eq. |6] The existence of such 
correlation is not surprising, since the radio and far-IR emission 
of star-forming galaxies are known to present a tight correlation 
(Yun et al. 2001, de Jong et al. 1985, Helou, Soifer & Rowan- 
Robinson 1985): the radio emission is predominantly produced 
by the synchrotron radiation of supernova remnants and the bulk 
of the far-IR emission is due to UV light from young and mas- 
sive stars reprocessed by interstellar dust. Hence, we consider 
this size estimate to be a good proxy for the global size of the 
star formation regions of galaxies. This is obviously an approx- 
imation, since this does not account for the dumpiness or gran- 
ularity of the region, but this is the best that we can do with 
existing datasets. 

Our local galaxy sample was cross-matched with the NRAO 
VLA Sky Survey (NVSS, Condon et al. 1998) and the Faint 
Images of the Radio Sky at Twenty-cm (FIRST, Becker, White 
& Helfand 1995), both obtained with the VLA at 20 cm. A total 
of 11, 47 and 58 galaxies have radio sizes in our ISO, AKARI 
and Spitzer local galaxy samples (see Table|2]). 

We computed the IR surface brightness using Eq. |7] 

Lir/2 

where the IR luminosity is divided by 2 since rjR is the far-IR 
(60 fim) half-light radius, which is derived from the 20 cm radio 
measurements using Eq. |6] 

5.1.2. Mid-IR compactness 

Using the low spectral resolution staring mode of the 
SpitzerflRS, Dfaz-Santos et al. (2010) measured the spatial ex- 
tent of the mid-IR continuum emission at 13.2yum for 211 lo- 
cal (U)LIRGs of the GOALS sample (see section 2.2.3). The 
13.2//m emission probes the warm dust (very small grains. 
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I„ = (L„/2)/(7iR„^) [Lo kpc-^] 



Fig. 11. Comparison of the two compactness indicators: Sir 
i-L}^/ (2X71 Rj^)), the IR surface brightness, and 13.2/zm com- 
pactness (percentage of unresolved Spitzer/IRS light at 13.2yL(m, 
Di'az-Santos et al. 2010). 



VSGs) heated by the UV continuum of young and massive stars, 
and hence traces regions of dust-obscured star formation. Instead 
of measuring the half-light radius of the sources at this wave- 
length, Diaz-Santos et al. (2010) calculated their fraction of ex- 
tended emission, or FEE, which they defined as the fraction of 
light in a galaxy that does not arise from its spatially unresolved 
central component. Conversely, the compactness of a source can 
be defined as the percentage of light that is unresolved, that is, 
lOOx(l-FEE). The angular resolution of Spitzer/IRS at 13.2/vm 
is ~ 3.6 " which, at the median distance of the sample used in 
this work, 91 Mpc, results in a spatial resolution of 1.7 kpc. 

In the following, we consider galaxies as "compact" if their 
13.2yum compactness is greater than 60%. With this definition, 
we find that 55% (1 17/21 1) of the GOALS galaxies are compact. 
Interestingly, while it is true that the fraction of galaxies show- 
ing compact star formation (i.e., compact hot dust emission) in- 
creases with increasing Lj°' (hence also with SFR), the compact 
population is not systematically associated with the most lumi- 
nous sources. On the contrary, galaxies with compact star forma- 
tion can be found at all luminosities (see Figure 4 of Dfaz-Santos 
et al. 2010). 

5.1 .3. Identification of tine galaxies with compact star 
formation 

In order to check whether both star formation compactness indi- 
cators are consistent, we used the 58 galaxies from the GOALS 
sample for which we can determine both Eir from the radio sizes 
(Sect. lSTOT l and a 13.2 yum compactness (Sect. l5?L2b . The com- 
parison of both compactness indicators shows a correlation with 
a dispersion of ~0.45 dex (Fig. fTTT l. The critical threshold of 
60 % in the 13.2 fim compactness above which we classify galax- 
ies as compact corresponds to Eir~3x10'" L© kpc"^. Hence, we 
hereafter classify as compact the galaxies for which Eir>3x10'" 
Lq kpc"^. This threshold is more than two orders of magnitude 
lower than typical upper limits for star formation on small (kpc) 
scales (see Soifer et al. 2001). We note that if it were not av- 
eraged on large scales, the local star formation surface density 
could be much higher in many of these sources. 
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Fig. 12. Distribution of IR sizes (half-light radius) of 119 lo- 
cal galaxies as derived from their 1.4 GHz radio continuum us- 
ing Eq. |6] Galaxies with an IR surface brightness greater than 
SiR=3xlO'" Lokpc"^, i.e., compact galaxies, are in blue, while 
extended galaxies are in red. The vertical dashed line indicates 
the typical resolution of ~ 1 .7 kpc of the mid-IR compactness at 
13.2jum estimated by Diaz-Santos et al. (2010). 

In Fig. [12l we present the distribution of far-IR sizes of ex- 
tended (red) and compact (blue) galaxies, estimated from radio 
20 cm imaging using Eq. |6] The median far-IR sizes of com- 
pact and extended galaxies are 0.5 kpc and 1.8 kpc respectively. 
The typical spatial resolution reached at 13.2jum, i.e., 1.7 kpc, 
is close to the typical size of extended galaxies and is signif- 
icantly larger than the median size for compact galaxies. This 
contributes to the relatively high dispersion seen in Fig. [TT| With 
a linear resolution of ~0.2 kpc at the average distance of the 
GOALS sample, the radio estimator is therefore a finer discrim- 
inant of compact galaxies when good quality radio data exist. 

In the following, we use both compactness indicators, i.e., 
radio and 13.2jum, with no distinction to define the projected 
IR surface brightness, Sir (=Lj°V(2x7r/?jj^)). For galaxies with 
a measured radio size, Rir is computed using Eq. |6] while for 
galaxies with a 13.2/im compactness estimate but no radio size 
we use the relation presented in Fig. [TT] 

5.2. IRS, a star formation compactness indicator 

The IRS ratio is compared to the IR surface brightness, Eir, in 
Fig-IE] The number of galaxies presented in this figure is larger 
than in Fig. [TT] because we include sources with no radio size 
estimate as well. We find that IRS is correlated with Zir for local 
galaxies following Eq. [8] 



0.15 



IRS = 0.22 [-0.05, -1-0.06] x 2^^' 



(8) 



where Sir is in Lq kpc"^. Hence IRS is a good proxy for the 
projected IR surface brightness of local galaxies. Galaxies with 
strong IRS ratios are also those which harbor the highest star- 
formation compactness. 

We showed in Fig. [TT]that galaxies having more than 60 % 
of their 13.2 fim emission unresolved by Spitzer-IRS, defined as 
compact star-forming galaxies by Diaz-Santos et al. (2010), pre- 
sented an IR surface brightness of Eir>3x10"' Lokpc"^. This 
threshold is illustrated in Fig. [13] by a vertical dotted line. It 
crosses the best-fitting relation of Eq. [8] at IRS-S, i.e., twice 
the central value of the Gaussian distribution of main sequence 
galaxies (Fig. [9] and Eq. HJ. As a result, compact star-forming 
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Fig. 13. Dependence of IRS (=l;°7L[8 //m]) with Sir 
(-L!^/(2x7tR^^)), the IR surface brightness. Galaxies to the 
right of the vertical dotted line are considered to be compact 
star-forming galaxies (Sir>3x10'" Lq kpc"^). Galaxies above 
the horizontal dotted line exhibit IRS ratio that is two times 
larger than that of main sequence galaxies, which follow the 
Gaussian IRS distribution shown in Fig.|9j. The sliding median 
is shown with open blue dots. It is fitted by the solid line in red 
and its 16th and 84th percentiles are fitted with the dashed red 
lines (see Eq.|8]i. 



galaxies, with Sir>3x10'" Lokpc"^, systematically present an 
excess in IRS, whereas nearly all galaxies with extended star 
formation exhibit a 'normal' IRS, i.e., within the Gaussian dis- 
tribution of Fig. [9] This is illustrated in Fig. [14] reproducing the 
IRS-LS diagram for local galaxies of Fig. [9] this time including 
the GOALS sample. The sub-sample of galaxies with measured 
IR surface brightnesses are represented with large symbols, with 
blue and red marking galaxies with extended and compact star- 
formation respectively, i.e.. Sir lower and greater than 3x10'" 
Lokpc"^. Galaxies with compact star formation systematically 
lie above the typical range of IRS values. The trend can be ex- 
tended to the local ULIRGs with no size measurement, since 
they are known to experience compact starbursts driven by ma- 
jor mergers (Armus et al. 1987, Sanders et al. 1988, Murphy et 
al. 1996, Veilleux, Kim & Sanders 2002). 

Very interestingly, the proportion of galaxies with compact 
star formation rises with Lj°' following a path very similar to 
the proportion of galaxies with IRS>S (Fig.[T5]l, the 68 % upper 
limit of the GOODS-//erac/ie/ galaxies (Eq.|3]l. Hence, IRS can 
be considered as a good proxy of the star formation compactness 
of local galaxies. This can be very useful for galaxies with no 
radio size measurement. 

In comparison, the fraction of excess IRS sources within 
the GOODS-//er5c/ie/ sample remains low (around 20%, see 
Fig. |6]-right and Fig. |9]l and never reaches such high propor- 
tions as seen in local ULIRGs. Due to the Herschel detection 
limit, however, only ULIRGs are individually detected at z>2. 
The IRS parameter is found to be biased towards high values in 
these galaxies which are responsible for the increase in the com- 
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Fig. 14. (top) and IR% ratio (bottom) versus L[8 /zm] for the 
local galaxy sample (ISO, AKARI, S'p/fzer-GOALS). Galaxies 
with compact mid-lR (more than 60% of the light emitted 
at 13.2 yum within the resolution element of 3.6") and radio 
(SiR>3xlO'" Lq) light distributions are marked with large filled 
red dots. Here the solid and dashed lines present the median and 
68 % dispersion around the Gaussian distribution as defined in 
EqH 
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Fig. 15. Fractions of compact local galaxies (black line), local 
galaxies with //?8>8 (blue dashed line) and GOODS-i/ersc/je/ 
galaxies with //?8>8 (red dash-dotted line) as a function of to- 
tal IR luminosity. The median redshift of the GOODS-Her^c/ie/ 
galaxies varies with L™, and is indicated by labels in grey along 
the red dash-dotted line. 



pactness fraction in the Herschel sample at the highest redshifts 
from 20 to 40%. 

Globally, this analysis suggests that compact sources have 
been a minor fraction of star-forming galaxies at all epochs, but 
locally, due to the low gas content of galaxies, compact sources 
make the dominant population of ULIRGs. Extending the anal- 
ysis of local galaxies to the distant ones, this also suggests that 
compact star formation takes place at all luminosities but does 
not dominate the majority of distant ULIRGs. Conversely, know- 
ing the compactness and mid-IR luminosity of a galaxy, one may 
optimize the determination of its total IR luminosity from mid- 
IR observations alone. This is discussed in Sect.|7] 

Finally, we have assumed in this section that the compact- 
ness measured either from the radio or from the mid-IR contin- 
uum is associated with star formation. We discuss the role of 
AGN in Sect. [8] but we can already note that when an AGN 
contributes to the IR emission of a galaxy, it does so mainly at 
wavelengths shorter than 20 ;um (Netzer et al. 2007, Mullaney et 
al. 2011a). If AGN were contributing to the infrared emission, 
they would tend to boost L% relative to Lj°', therefore reducing 
IR%. Instead, we see that an increasing compactness corresponds 
to an increase in //?8 as well. This reinforces the idea that we are 
dealing here with star formation compactness and not an effect 
produced by the presence of an active nucleus. In the next sec- 
tion, we show that compact star-forming galaxies are generally 
experiencing a starburst phase. 

5.3. IRE, a starburst indicator 

In the previous section, we have seen that high IR2, values were 
systematically found in galaxies with compact star formation re- 
gions. We now show that these galaxies are experiencing a star- 
burst phase. In the following, a star-forming galaxy is consid- 
ered to be experiencing a starburst phase if its "current SFR" is 
twice or more stronger than its "averaged past SFR" (<SFR>), 
i.e., if its birthrate parameter Z7=SFR/<SFR> (Kennicutt 1983) 
is greater than 2. Here <SFR>=M»/fgai, where fgai is the age of 
the galaxy. Alternatively, a star-forming galaxy may be defined 
as a starburst if the time it would take to produce its current stel- 
lar mass, hence its stellar mass doubling timescale, t, defined in 
Eq.lH 



T [Gyr] = [Mq] / SFR [Mq Gyr"'] = 1/sSFR [Gyr''] , (9) 

is small when compared to its age. Both definitions are equiva- 
lent if one assumes that galaxies at a given epoch have similar 
ages. 

In recent years, a tight correlation between SFR and Mt has 
been discovered which defines a typical specific SFR, (sSFR 
= SFR/M»), for "normal star-forming galaxies" as opposed to 
"starburst galaxies". This relation evolves with redshift but a 
tight correlation between SFR and M» is observed at all red- 
shifts from z~0 to 7 (Brinchmann et al. 2004, Noeske et al. 2007, 
Elbaz et al. 2007, Daddi et al. 2007a, 2009, Pannella et al. 2009, 
Magdis et al. 2010a, Gonzalez et al. 2011). Hence, we use the 
sSFR definition of a starburst since it can be applied at all look- 
back times. In the present section, we consider only local star- 
forming galaxies. The SFR - M, relation for local A/iTA/?/ galax- 
ies is shown in the left-hand part of Fig. [16] The best fit to this 
relation is a one-to-one correlation (0.26 dex-rms), hence a con- 
stant sSFR~0.25 Gyr ' or t~4 Gyr (Eq.|9]l. Local galaxies with 
compact star formation (large red dots), as defined in the previ- 
ous section, are systematically found to have higher sSFR, i.e., 
T< 1 Gyr, than that of normal star-forming galaxies. 
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Fig. 16. Left: SFR - M, correlation at z~0. Galaxies classified as compact are marked with large filled red dots. Solid line: fit to 
the main sequence SFR-M, relation: SFRocM*/[4xlO^ M©]. Dotted lines: 16th and 84th percentiles of the distribution around the 
sliding median (0.26 dex). Right:Relation of the sSFR and IR surface brightness of galaxies for which a radio size was estimated. 
The vertical dashed line illustrates the threshold above which galaxies have been classified as compact. The solid and dashed red 
lines are a fit to the sliding median of the relation (Eq. [TTt and its 68 % dispersion. 



Since there is a continuous distribution of galaxies ranging 
from the normal mode of star formation, with t~4 Gyr, to ex- 
treme starbursts, that can double their stellar masses in t~50 
Myr, we quantify the intensity of a starburst by the parameter 
RsB, which measures the excess in sSFR of a star-forming galaxy 
(which we label its "starburstiness"), as defined in Eq.fTOl 

RsB = sSFR/sSFRms = tms/t [> 2 for stai-bursts] , (10) 

where the subscript MS indicates the typical value for main 
sequence galaxies at the redshift of the galaxy in question. A 
starburst is defined to be a galaxy with Rsb>2. 75 % of the 
galaxies with compact star formation (Sir>3x10'" Lq kpc"^) 
have /?SB>2, hence are also in a starburst mode, and 93% of 
them have /?sb>1- Conversely, 79% of the starburst galaxies 
are "compact". Globally, starburst galaxies with sSFR > 2 x 
<sSFR> have a median Sir~1.6x10^' L© kpc"^, hence more 
than 5 times higher than the critical IR surface brightness above 
which galaxies are compact. The size of their star-forming re- 
gions is typically 2.3 times smaller than that of galaxies with 
sSFRms. 

The sSFR and Sir are correlated with a 0.2 dex dispersion 
(Fig. [T6l-right) following Eq. [TT] where sSFR is in Gyr ' and 

IiR in Lokpc^^: 

sSFR= 1.81 [-0.66, +1.05] X lO"'* x I.^^^ (11) 

Both parameters measure specific quantities related to the SFR: 
the sSFR is measured per unit stellar mass, while Sir is related 
to the SFR (derived from Lj°') per unit area. Note, however, that 
it was not obvious a priori that these quantities should be corre- 
lated, since the stellar mass of most galaxies is dominated by the 
old stellar population, whereas the IR (or radio) size used to de- 
rive Sir measures the spatial distribution of young and massive 
stars. 

Because the starburstiness and the /7?8 ratio are both en- 
hanced in compact star-forming galaxies, they are also correlated 
as shown by Fig.[T7l The fit to this correlation is given in Eq.fTSl 

RsB = (//?8/4)'-2 (12) 



lOOF 




Fig. 17. /;sb=sSFR/sSFRms versus IR8 {^H^/Lf.) forz~0 galax- 
ies {AKARI and GOALS samples). The red line is the best fit 
(plain) and its 0.3 dex dispersion (dashed). 



The dispersion in this relation is 0.3 dex. Hence we find that it 
is mostly compact starbursting galaxies that present atypically 
strong //?8 bolometric correction factors, although there is not 
a sharp separation of both regimes, but instead a continuum of 
values. 



6. IRS as a tracer of star formation compactness 
and "starburstiness" in distant galaxies 

In the previous section, we have defined two modes of star for- 
mation: 

- a normal mode that we called the infrared main sequence, in 
which galaxies present a universal IRE bolometric correction 
factor and a moderate star formation compactness, Sir, and 
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Fig. 18. Redshlft evolution of the median specific SFR 
(sSFR=SFR/M,) of star-forming galaxies. Values for individual 
GOODS-i/eric/ie/ galaxies are shown as grey points. Median 
sSFR values in redshift bins are shown with open circles (blue 
for GOODS-N and black for GOODS-S). Values combin- 
ing individual detections and stacking measurements for unde- 
tected sources are shown with filled triangles (blue upward for 
GOODS-N and black downward for GOODS-S). The red solid 
line is the fit shown in Eq. [13] and the dashed lines are a fac- 
tor 2 above and below this fit. Starbursts are defined as galaxies 
with a sSFR>2xsSFRms (blue zone). The yellow zone shows the 
galaxies with significantly lower sSFR values. 

- a starburst mode, identified by an excess SFR per unit stellar 
mass, hence sSFR, as compared to the typical sSFR of most 
local galaxies. 

Galaxies with an enhanced IR2i ratio were systematically found 
to be forming their stars in the starburst mode and to show a 
strong star formation compactness (Sir>3x10"* Lokpc"^). In 
order to separate these two modes of star formation in distant 
galaxies as well, we first need to define the typical sSFR of 
star-forming galaxies in a given redshift domain. This defini- 
tion has become possible since the recent discovery that star- 
forming galaxies follow a tight correlation between their SFR 
and M» with a typical dispersion of 0.3 dex over a large range 
of redshifts: z~0 (Brinchmann et al. 2004), z~l (Noeske et al. 
2007, Elbaz et al. 2007), z~2 (Daddi et al. 2007a, Pannella et al. 
2009), z~2> (Magdis et al. 2010a), z~4 (Daddi et al. 2009, Lee et 
al. 201 1) and even up to z~l (Gonzalez et al. 2011). 

6.1. Evolution of the specific SFR witfi cosmic time and 
definition of main sequence versus starburst galaxies 

In the following, we assume that the slope of the SFR - M» re- 
lation is equal to 1 at all redshifts, hence that the specific SFR, 
sSFR (=SFR/M,), is independent of stellar mass at fixed red- 
shift. A small departure from this value would not strongly affect 
our conclusions and the same logic may be applied for a differ- 
ent slope. At z~0, our local reference sample is well fitted by a 
constant sSFR (see Fig. [T6l-left). although the best-fitting slope 
is 0.77 (Elbaz et al. 2007). At z~l+0.3, Elbaz et al. (2007) find 
a slope of 0.9 but we checked that the dispersion of the data al- 
lows a nearly equally good fit with a slope of 1 . At z~2, Pannella 
et al. (2009) find a slope of 0.95 consistent with the value ob- 
tained by Daddi et al. (2007a) in the same redshift range. Lyman- 



break galaxies at z~3 (Magdis et al. 2010a) and z~4 (Daddi et 
al. 2009) are also consistent with a slope of unity. From a differ- 
ent perspective, Peng et al. (2010a) argue that a slope of unity 
is required to keep an invariant Schechter function for the stellar 
mass function of star-forming galaxies from z~0 to 1 as observed 
from COSMOS data, while non-zero values would result in a 
change of the faint-end slope of the mass function that would be 
inconsistent with the observations. 

However, the slope of the SFR - M, relation is sensitive to 
the technique used to select the sample of star-forming galaxies. 
Karim et al. (2011) find two different slopes depending on the 
selection of their sample: a slope lower than 1 for a mildly star- 
forming sample, and a slope of unity when selecting more ac- 
tively star-forming galaxies (see their Fig. 13). Using shallower 
Herschel data than the present observations, Rodighiero et al. 
(2010) found a slope lower than unity. 

Assuming that the slope of the SFR - M, relation remains 
equal to 1 at all redshifts, a main sequence mode of star forma- 
tion can be defined by the median sSFR in a given redshift in- 
terval, sSFRms(z). The starburstiness, described in Eq. [10] mea- 
sures the offset relative to this typical sSFR. Since at any redshift 
- at least in the redshift range of interest here, i.e., z<3 - most 
galaxies belong to the main sequence in SFR - M,, we assume 
that the median sSFR measured within a given redshift interval 
is a good proxy to the sSFRms(z) defining the MS. Galaxies de- 
tected with Herschel follow the trend shown with open circles in 
Fig.[T8](blue for GOODS-N and black for GOODS-S). We have 
performed the analysis independently for both GOODS fields 
in order to check the impact of cosmic variance on our result. 
To correct for incompleteness, we performed stacking measure- 
ments as for Fig. [6] but in redshift intervals. The stacking was 
done on the PACS-lOOyum images using the 24/im sources as a 
list of prior positions. The resulting values (blue upward trian- 
gles for GOODS-N and black downward triangles for GOODS- 
S) were computed by weighting detections and stacking mea- 
surements by the number of sources used in both samples per 
redshift interval. The SFR was derived from Lj°' extrapolated 
from the PACS-100//m band photometry using the CEOl tech- 
nique. The CEOl method works well for 100 /vm measurements 
up to z~3 as noted already in Elbaz et al. (2010), and we confirm 
this agreement with the extended sample of detected sources 
in the present analysis (Sect. I7.3l l. The trends found for both 
fields are in good agreement. The stacking + detection measure- 
ments for GOODS-N are slightly lower than those obtained for 
GOODS-S which may result from a combination of cosmic vari- 
ance and the fact that the GOODS-S image is deeper 

The redshift evolution of sSFRms(z) (Fig. [TSl ). accounting 
for both detections and stacked measurements, is well fitted by 
Eq.[l3] 

sSFRu^ [Gyr-'] = 26 x t'l^^^ , (13) 

where fcosmic is the cosmic time elapsed since the Big Bang in 
Gyr. A starburst can be defined by its sSFR following Eq. [14] 

sS FRsB [Gyr-'] > 52 xCL. (14) 

The intensity of such starbursts, or "starburstiness", is then de- 
fined by the excess sSFR: /?sb=sSFRsb/sSFRms. Due to the evo- 
lution observed with cosmic time, a galaxy with a sSFR twice as 
large as the local MS value would be considered a starburst to- 
day, but a galaxy with the same sSFR at z~ 1 would be part of the 
main sequence. 

We have seen that for local galaxies, the starburstiness and 
IR2, are correlated (see Fig. [TTt . The same exercise for distant 
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Fig. 19. T^sb^sSFR/sSFRms versus im (=L}^/U) for the distant 
GOODS-He rschel galaxies. The green lines show the range of 
values occupied by main sequence galaxies (68 % dispersion) 
in sSFR (horizontal lines; /?sb = 1±1) and //?8 (vertical lines; 
//?8ms=4±1.6, see Eq. |5]). The thick grey line show the sliding 
median for the GOODS-Herschel galaxies. The diagonal blue 
lines are the best fit (solid) and 68 % dispersion (dashed) for lo- 
cal galaxies as in Fig. [17] Large open dots show the position of 
sub-mm galaxies from Menendez-Delmestre et al. (2009) and 
Pope et al. (2008a). 
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GOODS-Herschel galaxies, mixing galaxies of all luminosities 
and redshifts, is shown in Fig. [19] Distant galaxies exhibit a non 
negligible dispersion, but their sliding median, shown by a thick 
grey line in Fig. [19] is coincident with the best fit relation for 
local galaxies (solid and dashed blue lines). 

We find that 80 % of the galaxies which belong to the SFR 
- M» main sequence - with 0.5<^sb^2 - also belong to the 
main sequence in IR2i - //?8ms=4+1.6 (Eq. O. Hence we con- 
firm that the two definitions of "main sequence galaxies" are 
similar and that on average they represent the same galaxy pop- 
ulation. We note also that even though there is a tail toward 
stronger starburstiness and compactness, i.e., increased 7?sb and 
//?8, this regime of parameter space is only sparsely populated 
in the GOODS-//er5c/ie/ sample, which suggests that analogs 
to the local compact starbursts predominantly produced by ma- 
jor mergers remain a minority among the distant galaxy popula- 
tion. Finally, we see that sub-mm galaxies (large open circles in 
Fig. [19]) also follow the same trend. 

6.2. Star formation compactness of distant galaxies 

We have shown that local galaxies with high Sir values also ex- 
hibit high IRdi ratios. We do not have IR or radio size estimates 
for the distant galaxy population, but we can use the high reso- 
lution HST-ACS images to study the spatial distribution of the 
rest-frame UV light in the populations of MS and SB galaxies. It 
has been suggested that distant (U)LlRGs at 1.5<z<2.5 (Daddi 
et al. 2007a) and at z~3 (Magdis et al. 2010b) are not optically 
thick since the SFR derived from the UV after correcting for ex- 
tinction using the Calzetti et al. (2000) law is consistent with the 
SFR derived from radio stacking measurements at these redshifts 
(see also Nordon et al. 2010). 



Fig. 20. Stacked images (5" on a side) centered on main se- 
quence (left column) and starburst (right column) galaxy posi- 
tions. Typical MS galaxies are selected to have 7?sb = 1±0.1 and 
SB galaxies /?sb^3. Each image results from the stacking of 
HST-ACS images in the B (4350 A), V (6060 A) and 1 (7750 A) 
bands corresponding to the rest-frame UV at ~2700A at z=0.7 
(first line), z=1.2 (second line) and z=1.8 (third line) respec- 
tively. 

Table 3. UV - 2700 A half-light radii of distant main sequence 
and starburst galaxies. 



Redshift 


Main Sequence 


Starburst 






^SB>2 ^SB>3 


0.7 


5.2 kpc 


3.9 kpc 2.5 kpc 


1.2 


4.4 kpc 


3.3 kpc 2.5 kpc 


1.8 


3.0 kpc 


2.5 kpc 2.0 kpc 



We use HST-ACS images in the B (4350 A), V (6060 A) and 
/ (7750 A) bands to sample the same rest-frame UV wavelength 
of ~2700A at z=0.7, 1.2 and 1.8 respectively. MS galaxies are 
selected to have 7?sb = 1±0.1 (Eq. [TOli whereas SB galaxies are 
defined as galaxies with 7?sb^2. We also tested a stricter def- 
inition for starbursts, 7?sb^3 to avoid contamination from MS 
galaxies (Table O. The result of the stacking of HST-ACS sub- 
images is shown in Fig.|20]for MS (left column) and SB galaxies 
with /?SB^3 (right column). It is clear that the sizes of the star- 
bursts are more compact than those of the main sequence galax- 
ies. The half-light radius of each stacked image was measured 
with GALFIT (Peng et al. 2010b) and is Hsted in Table[3] 

These sizes are consistent with those obtained by Ferguson et 
al. (2004). SB galaxies typically exhibit half-light radii that are 
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two times smaller than those of MS galaxies, implying projected 
star formation densities that are 4 times larger We verified that 
this was not due to a mass selection effect by matching the stel- 
lar masses in both samples and obtained similar results, although 
with larger uncertainties. These sizes are larger than the radio- 
derived IR half-light radii of the local sample of MS (1.8 kpc) 
and SB (0.5 kpc) galaxies. However, since the distant galaxy 
sample has a different mass and luminosity selection than that 
of the local reference sample, we cannot directly compare their 
sizes. However, the difference in the relative sizes among the 
high-redshift galaxies confirms that star formation in distant star- 
bursts is more concentrated than that in distant main sequence 
galaxies. This, again, is strong evidence for a greater concen- 
tration of star formation in galaxies with higher specific SFRs. 
Since we have seen that sSFR and //?8 are correlated (Fig. [19}, 
this implies that in distant galaxies, like in local ones, galaxies 
with strong //?8 ratios are likely to be compact starbursts. 

This result is consistent with the work of Rujopakam et al. 
(201 1), who measured IR luminosity surface densities for distant 
(U)LIRGs similar to those found in local normal star-forming 
galaxies. However, we find that this is the case for most but not 
all high redshift (U)LIRGs. Compact starbursts do exist in the 
distant Universe, even among (U)LlRGs, but they are not the 
dominant population. 

7. Toward a universal IR SED for Main Sequence 
and Starburst galaxies 

7.1. Medium resolution IR SED for main sequence (IR8~ 4) 
and starburst (IRS> 8j galaxies 

At z<2.5 - where we can estimate the rest-frame Lg from 
Spitzer IRAC, IRS and MIPS photometry as well as reliable 
LiR from Herschel measurements at rest-frame /l>30jum - the 
IR% i-LYnlL'i) ratio follows a Gaussian distribution centered on 
//?8~4 (Eq.|5j Fig.|9]l, with a tail skewed toward higher values for 
compact starbursts. This defines two populations of star-forming 
galaxies or, more precisely, two modes of star formation: the 
MS and SB modes. Galaxies in the MS mode form the Gaussian 
part of the IR% distribution and present typical sSFR values (i.e., 
^SB~1) while SB exhibit stronger IR% values (see Fig.|9j and a 
stronger "starburstiness" (/?sb>2). 

/7?8 is universal among MS galaxies of all luminosities and 
redshifts. This suggests that these galaxies share a common IR 
SED. In the local Universe, the rest-frame L\2, L25, Leo, ^loo 
from IRAS and L15 from ISOCAM were also found to be nearly 
directly proportional to Lj"' (see CEOl and Elbaz et al. 2002), 
hence reinforcing this idea. To produce the typical IR SED of 
MS and SB galaxies, we use ^-correction as a spectroscopic tool. 
We separate MS and SB galaxies by their IR8 ratios: IR8-4+2 
for MS galaxies (as in Eq.jSj and /7?8>8 (hence >2cr away from 
the MS) for SB galaxies. We then normalize the individual IR 
SEDs by a factor 10' '/^ir so that all galaxies are normalized 
to the same reference luminosity of Lj°'=10" L©. The result is 
shown with light grey dots in the left-hand part of Fig. |2T]for MS 
galaxies and in the right-hand part of Fig. |2T]for SB galaxies. 
A sliding median was computed in wavelength intervals which 
always encompass 25+5 galaxies (blue points for MS in Fig.|2T} 
left and red points for SB in Fig.|2T}right). As a result, the typical 
MS and SB IR SEDs have an effective resolution of A/AA=25 and 
10 respectively, nearly homogeneously distributed in wavelength 
from 3 to 350 /im. 

The typical MS IR SED in the left-hand part of Fig. [21] has 
a broad far-IR bump centered around 90 fim, suggesting a wide 



range of dust temperatures around an effective value of -30 K, 
and strong PAH features in emission. Instead, the typical IR SED 
for SB galaxies (Fig.|2T}right) presents a narrower far-IR bump 
peaking around /l~70-80 fim, corresponding to an effective dust 
temperature of -40 K, and weak PAH emission lines. We note 
however, that these prototypical IR SEDs result from the com- 
bination of 267 and 111 galaxies for the MS and SB modes, 
respectively. They therefore should be considered as average 
SEDs, acknowledging that there is a continuous transition from 
one to the other with increasing IRE or star-formation compact- 
ness. In the next Section, we provide a model fit to these SEDs 
to better describe their properties. 

7.2. SED decomposition of main sequence and starburst 
galaxies 

In order to interpret the physical nature of the MS and SB SEDs 
derived in the previous section, we adopt a simple phenomeno- 
logical approach. We decompose the two classes of SEDs with 
the linear combination of two templates, shown in Fig.|22l (1) 
a ''star-forming region" component including Hii regions and 
the surrounding photo-dissociation region (labeled SF), and (2) 
a "diffuse ISM" (interstellar medium) component accounting for 
the quiescent regions (labeled ISM). The luminosity ratio of the 
two components controls the IRS parameter This SED decom- 
position is not unique and the two components used here are not 
rigorously associated with physical regions of the galaxies. 

The SED of each sub-component is given by the model of 
Galliano et al. (201 1, in prep.; also presented by Galametz et al. 
2009). This model adopts the Galactic dust properties of Zubko, 
Dwek & Arendt (2004). To account for the diversity of physical 
conditions within a galaxy, we combine the emission of grains 
exposed to different starlight intensities, U (normalized to the 
solar neighborhood value of 2.2 x 10"^ Wm"^). We assume, 
following Dale et al. (2001), that the mass fraction of dust ex- 
posed to a given starlight intensity follows a power-law (index 
Of): dMdust/dU oc U^". The two cutoffs are and t/,„,„ + At/. 
We fit the two SEDs simultaneously, varying only the luminos- 
ity ratio of the two components. We add a stellar continuum to 
fit the short wavelengths (see Galametz et al. 2009 for a descrip- 
tion). This component is a minor correction. In summary, the 
free parameters for the fit are: 

- the starlight intensity distribution parameters (a, Umin and 
AU) of each sub-component; 

- the PAH mass fraction and charge of each sub-component; 

- the luminosity ratio of the two components for the main se- 
quence and for the starburst; 

- the contribution of the stellar continuum (negligible here) 

The fits are shown with solid black lines in Fig.|2T]while the 
derived templates that we used for the decomposition are shown 
in blue and red lines in Fig. [22] The most relevant parameters are 
summarized in Table|4| The "diffuse ISAf' component has colder 
dust and a larger PAH mass fraction than the "star-forming re- 
gion" SED. 

The main differences between galaxies in the MS and SB 
modes are: 

- the effective Tjust of galaxies in the SB mode is warmer than 
that of MS galaxies, i.e., ~40 K versus ~31 K; 

- the contribution of diffuse ISM emission to the SB SED 
is negligible (8 %), consistent with the strong compactness 
seen both for local starbursts (in radio and mid-IR imaging. 
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Fig. 21. Composite spectral energy distribution of the typical main sequence galaxy {left; IRdi-4+2, see Eq.|5]) and starburst (right; 
//?8>8, i.e., above 2cr). Light grey dots: individual GOODS-Herschel galaxies normalized to Lj°* = 10" Lq. The large filled 
symbols with error bars are the median and associated uncertainty of the MS (left figure, blue dots) and SB (right figure, red dots) 
galaxies computed in intervals of wavelengths defined to contain a fixed number of 25+5 galaxies. The uncertainty on the median 
values is derived from the 16th and 84th percentiles around the median divided by the square root of the number of galaxies. The 
model fit to each SED is shown with a solid black line while the opposing SED (MS or SB) is shown with a dotted black Une for 
comparison. 
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Fig. 22. Components used in the fit ofFig.\2l\ These two com- 
ponents have been constrained by the simultaneous fit of the two 
SEDs (Fig.l2TI). The main sequence and starburst SEDs are the 
linear composition of these two components. The luminosity ra- 
tio of these components controls IRS. 



see Sect. |5]l and for high redshift analogs (in the rest-frame 
UV, see Sect.lOli: 
- the MS SED requires a wider distribution of dust tempera- 
tures, typically ranging from 15 to 50 K; 



- the stronger contribution of PAH lines to the broadband mid- 
IR emission in the MS SED is the main cause for the differ- 
ence in IR2t ratios between the two populations. 

Galaxies are distributed continuously between the MS and vari- 
ous degrees of SB strength, hence this decomposition technique 
can be used in the future to produce SEDs suitable for ranges 
of IRE or sSFR values, in the form of a new library of tem- 
plate SEDs. We note, however, that this decomposition of the 
typical MS and SB SEDs is not unique. For example, the SB 
SED is very similar to the CEOl template for a local galaxy with 
LiR = 6 X 10" Lq galaxy in the local Universe, which turns out 
to be close to the observed median luminosity of the starbursts. 
Instead, the MS SED is closer to the CEOl SED for a 4x10'^ Lq 
galaxy in the local Universe. 

We note also that a direct fit of the Rayleigh- Jeans portions 
of both SEDs would favor an effective emissivity index of jS=1.5 
for the MS and /3-2 for the SB. However, this is a degenerate 
problem. Indeed, the effective emissivity index /3 is not neces- 
sarily equal to the intrinsic /3 of the grains. A temperature dis- 
tribution of grains having an intrinsic fi - 2 would flatten the 
sub-mm SED and can give an effective /3 of ^ 1.5, as it is the 
case for our star-forming region. Finally, it is also not possible 
to disentangle some potential contribution from an AGN, par- 
ticularly for the SB SED. Indeed, AGN are known to be ubiq- 
uitous in LIRGs (Iwasawa et al. 201 1) and ULIRGs (Nardini et 
al. 2010), and they may contribute in part to the mid-IR con- 
tinuum, mostly in SB SEDs, since those are both more compact 
and exhibit lower PAH equivalent widths than they do MS galax- 
ies. However, even if AGN may contribute to some fraction of 
the light in these galaxies, they cannot dominate both in the mid 
and far-IR regimes since we find evidence that PAHs dominate 
around 8 fim in both MS and SB galaxy types, even if they are 
stronger in the MS SED. The high IR2, values measured for SBs 
also suggest that star formation dominates the IR emission in 
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Table 4. Main parameters relative to the SED decomposition. 



Dust Conditions 




ISM 


Star Forming 


[units] 
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[2.2 X IQ-' Wm-^] 




19 


53 


[K] 


/pah 
IRS 


1.9 
15 


0.3 
3 


[4.6 %] 


SED properties 




Main Sequence 


Starburst 




'P 


62 


92 


[%] 


IRS 


5 


11 




rppcak 


31 


40 


[K] 



10 r 




Notes. (U) is the luminosity averaged starlight intensity; 
Tcs - (f/)''* X 17.5 K is the corresponding effective tempera- 
ture of the grains; /pah is the PAH-to-total-dust mass fraction; 
(/> = Lsf/{Lkm + ^sf) is the luminosity fraction of the star-forming 



component. T'^^'^^ is the effective dust temperature corresponding to the 
peak of the far-lR bump using Wien's law. 
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these galaxies. In Sect. [8] we present a technique to search for 
hidden AGN activity in the GOODS-//er.sc/ie/ galaxies. 



7.3. Derivation of total IR luminosities from monochromatic 
measurements 

Now that we have defined a typical IR SED for main sequence 
galaxies, this SED may be used to extrapolate the total IR lu- 
minosity of galaxies for which only one measurement exists. 
Ideally, one would need to know the value of IRE or equiva- 
lently the starburstiness, Rs,b, of a galaxy, to know whether to 
use the MS or SB SED. But this would require to already know 
the actual SFR of a galaxy, which is what we are looking for An 
alternative technique would consist in using the star formation 
compactness of a galaxy, or Eir, to determine if it is in the main 
sequence or starburst mode. 

Assuming that all galaxies share the same MS SED, in 
Fig. |23]we compare the total IR luminosity that can be extrap- 
olated from a single passband using the main sequence IR SED 
(Lj'j^) with the value (L!^) measured from an SED fit to galax- 
ies with "clean" Herschel detections in several bandpasses (see 
Sect. 12.31 1. Both luminosities agree with an average uncertainty 
of ~35 % when L^^ is extrapolated from 24 yum and -20 % when 
L^j^ is derived from one of the 100 to 350 ^m. wavelengths. This 
is remarkable since we only used di single IR SED to extrapo- since Tfjgg- 

fall 1/ y[Gp 

gas+stars 

late L^^ for all galaxies. The MS SED does a better job than the 
CEOl technique (see Fig. 3 in Elbaz et al. 2010), which overes- 
timates Lj°' from lAjjm measurements at z>1.5 as well as from 
SPIRE 250 and 350 jum measurements at z<1.3. Note however 
that individual galaxies do present a wide range of IR SEDs with 
different dust temperatures. 

Finally, we note that these extrapolations work nearly 
equally well for X-ray AGN (black open triangles in Fig. |23) on 
average, although the dispersion is slightly larger for these galax- 
ies. This suggests that star formation dominates the IR emission 
in the hosts of typical AGN in deep-field X-ray surveys. We dis- 
cuss the properties of AGN in detail in Sect. |8] 




1 1.5 

redshift 



2.5 



Fig. 23. Ratio of the extrapolated {L'^^ over Herschel total IR 
luminosity as a function of redshift, for all clean GOODS- 
Herschel galaxies. L^^ is computed by normalizing the main 
sequence SED to the broadband photometric measurement at 
A. The 5 passbands used for the extrapolation are, from the 
top to bottom, Spitzer MIPS-24//m, Herschel PACS-lOOjum 
& 160 i^m, Herschel SPIRE-250yum & 350yum. LiR^Lf^'^''^'^*"^' is 
measured using the full set of Herschel measurements at rest- 
frame wavelengths /l>30;um to normalize the main sequence 
SED and by integrating over 8-1000 ;um. Black triangles: AGN. 
The solid lines with error bars are the sliding median and the 
16th and 84th percentiles around it. Upper panel: SED of M82 
and MIPS 24yL(m filter at z=0.25, 0.9, 1.4, 2. 



7.4. Interpretation of the connection between compactness, 
starburstiness and IRE 

At fixed redshift, a normal galaxy forms stars at a rate propor- 
tional to its gas mass divided by the free-fall time, as expressed 
inEq.ini 



SFR oc Mgas/Tfree-fall Mg^s Pg^^s+stais 



(15) 



If one assumes that the free-fall time is dominated by the gas 
density, then it follows that the right term of the equation is pro- 
portional to Pgfs, as in the Schmidt- Kennicutt relation and close 
to the value of 1.4 found by Kennicutt (1998b) for projected gas 
and SFR densities. However the role of stars (in the free-fall 
time) may not be negligible in some conditions and may partly 
explain why including them in the relation may reduce the ob- 
served dispersion, as proposed by Shi et al. (201 1). 

If, instead, we consider separately the roles played by the 
density and gas mass in Eq. [15] we may interpret that SB galax- 
ies form stars more efficiently as a result of a greater Pgas+stars, 
hence shorter Tfree-faii, possibly due to a merger In the case of 
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more distant galaxies, a galaxy with a similar stellar mass will 
naturally possess a higher gas fraction, and hence gas mass for 
its stellar mass, which will result in a greater sSFR. In this frame- 
work, where the difference of sSFR with redshift comes from a 
greater gas mass content in the past, but with similar gas den- 
sities, it is natural that MS galaxies exhibit similar IRE values 
and share a common prototypical IR SED. In the case of a SB, 
where the density is increased (e.g., by a merger), the IRE ratio 
is very sensitive to the geometry of the young stellar population 
(see, e.g., Galliano, Dwek & Chanial 2008, in particular their 
Fig. 6). Increasing the compactness of the young stellar popula- 
tion, hence Eir, in the case of SB would increase the radiation 
field and push the photo-dissociation region farther away where 
molecules such as PAHs can survive and emit their light. The 
equivalent width of PAHs would then be reduced and the contri- 
bution of continuum emission increased, resulting in greater IRE 
values. 

8. Unveiling dusty AGN within starburst galaxies 

8.1. The IRE bolometric correction factor for X-ray and 
power-law AGN 

It is known that AGN can heat the dust that surrounds them 
to temperatures of several hundreds of Kelvin and produce a 
spectrum that can dominate the mid-IR emission of a galaxy, 
but which typically falls off steeply at wavelengths longer than 
20yum (Netzer et al. 2007; see also MuUaney et al. 201 la). AGN 
can therefore make quite a significant contribution to the emis- 
sion around 8 //m. In this study, we have identified AGN using all 
available criteria and have excluded them from most aspects of 
the analysis up to this point. Only highly obscured and unidenti- 
fied AGN may still be present in the sample that we used to carry 
out the IRE analysis. 

Before trying to define techniques to identify those hidden 
AGN, we should examine the properties of known and well- 
recognized AGN in the GOODS-Herschel data. We divide the 
known AGN into two populations (defined in Sect. 13. Il l: the X- 
ray/optical AGN and the infrared power-law AGN. The reason 
for this separation is that infrared power-law sources already 
show evidence that part of their mid-IR emission is powered by 
an AGN, since they have been identified as galaxies with a rising 
mid-infrared continuum from 3.6 to 8 /vm. We used the criteria 
given in Eq.[T6l taken from Le Floc'h et al. (in prep.; technique 
similar to that from Ivison et al. 2004 and Pope et al. 2008b for 
sub-mm galaxies): 

Sy[4.5nm] > Sv[3.6fj.m] 

'«^>o (^)< 3.64 xZ.,,o(£f^)-. 0.15 (16) 

X-ray/optical AGN that also meet the infrared power-law criteria 
are counted as power-law sources in the following discussion. 

Surprisingly, AGN of both types exhibit IRE ratios that are 
centered on the median of star-forming galaxies (upper panels of 
Fig.l24li. i.e., IRE~4.9 (Eq.O. This behavior may be understood 
for the non-power-law X-ray/optical AGN, since there is no ev- 
idence that high amounts of radiation from the AGN heats the 
surrounding dust in these galaxies. Hence for those galaxies, it 
is the star formation that is most probably responsible for both 
the 8 jum and far-IR emission. 

The situation is different for power-law AGN, however In 
this case, we know, by definition, that a hot dust continuum is 
present that exceeds the stellar continuum emission at wave- 
lengths shorter than 8 //m. Still, the IRE ratios observed for these 



galaxies remains similar to that found for star-forming galaxies 
(upper panel of Fig. l24l l. This is consistent with Fig. |23] where 
we showed that extrapolations of the total IR luminosity from 
any single photometric measurement between the observed 24 
to 350yum passbands were nearly as accurate for AGN as for 
star-forming galaxies with no X-ray or optical AGN signatures. 

However, we have seen in Sect. |6] that compact starbursts 
have larger IRE ratios than do normal star-forming galaxies. 
Hence it is possible that two mechanisms act in opposite ways: 
some contribution from the hot dust heated by an AGN to the 
mid-IR light may be counterbalanced by the presence of a star- 
burst that increases the far-IR over mid-IR ratio. To test this pos- 
sibility, we correct for the effect of starbursts on IRE in the next 
section. 



8.2. Correction of tine effect of starbursts on IRE 

We showed that a starburst induces an enhancement of the far-IR 
emission at fixed 8 fim luminosity, whereas an AGN may induce 
an increase of the Ej-im luminosity at fixed far-IR luminosity. 
The enhancement of IRE in the presence of a starburst is propor- 
tional to its intensity, as measured by the starburstiness, 7?sb (see 
Fig. [TTland Fig. [T9] l. Hence it can be corrected by normalizing 
IRE by RsB, i e., replacing IRE by IRE/Rsa. 

The de-boosted IRE ratios galaxies with known AGN 
are shown in the lower panels of Fig. |24] (open triangles). 
Interestingly, we find that the two AGN populations behave dif- 
ferently. The IRE ratios of X-ray/optical AGN (from which we 
have excluded power-law AGN) remain centered on the region 
of star-forming galaxies. The fraction of AGN falling below the 
lower limit of main sequence star-forming galaxies increases 
from 1 1 % to 22 % but the same happens above the upper limit, 
which illustrates that this is just a result of the enhanced disper- 
sion produced when correcting by /?sb- This suggests that the 
IR emission of X-ray AGN is predominantly powered by star 
formation (as also confirmed by Mullaney et al. 2001b). 

The case of the power-law AGN is very different. The frac- 
tion of galaxies falling below the lower limit in IRE of MS star- 
forming galaxies rises from 33 % (already three times higher 
than for the X-ray AGN) to 70 % after dividing by /?sb- Hence 
the majority of the power-law AGN show evidence for an 8 fim 
excess, but this excess mid-infrared emission was disguised by 
the presence of a concurrent starburst. 

Two important conclusions can be derived from this obser- 
vation: 

1 . most of the IR emission from non-power-law X-ray/optical 
AGN appears to be powered by dust heated by stars. Hence 
their IR luminosities may be used to derive SFRs; 

2. the bulk of power-law AGN host both an obscured AGN and 
a compact starburst. 

The second point suggests that there is a physical link be- 
tween both activities, the obscured AGN and the starburst, since 
they take place at the same time. It makes sense that infrared 
power-law AGN are associated with both compact starbursts and 
obscured active nuclei since compactness is required both to ex- 
plain the excess sSFR of these galaxies as well as their dust 
obscuration. We note that the power-law criterion has not been 
demonstrated to be a perfect tracer of dusty AGNs, hence some 
of the galaxies selected by this criterion may just be purely star- 
forming galaxies. 



D. Elbaz et al.: GOODS-Herschel: an infrared main sequence for star-forming galaxies 



23 



100 



_T 10 



0.1 
100 



10 



CD 

cr 1 



0.1 



A tA, - A- 



0.5 1.5 2.5 



-H+l 1 1 — I I I I I I I 1 1 — I I I I I I I 

V V V 



V ▼ 

V- 

V y TO. 



V 

V V 



V * ▼ 



10^ 



10^ 



100 



J* 10 



00 



0.1 
100 



10 



oo 

q; 1 



0.1 



"T^a" 



vr 



A A 

Z 



0.5 1.5 2.5 



H — I I I I I I 



V 



10^ 



10^ 



10' 



La (Lo) 



Fig. 24. Left: IR?i=L}^/Lfi, ratio as a function of Lg for X-ray/optical AGN excluding power-law AGN. Upper panel: As-measured 
/7?8 ratios of the X-ray/optical AGN. They fall on the same sequence as star-forming galaxies. Bottom panel: The same galaxies 
after de-boosting their //?8 ratios for starburstiness. The galaxies present a slightly different configuration, but basically remain 
centered on the trend of star-forming galaxies, suggesting that in these AGN the IR emission is dominated by star formation. Right: 
IR% ratio as a function of Lg for power-law AGN. Upper panel: As-measured IR% ratios of the power-law AGN. They fall on the 
same sequence as star-forming galaxies, like the X-ray /optical AGN at left. Bottom panel: The same galaxies after de-boosting 
their IR% ratios for starburstiness. A significant fraction of the galaxies now fall systematically below the main sequence IR2> ratio, 
showing that their 8 pm emission is boosted by the hot dust heated by the AGN. 



8.3. Searching for unknown obscured AGN 

We have seen in Sect. I8.2l that the presence of a starburst could 
hide the signature of a dusty AGN on //?8, and that this could be 
corrected by normalizing //?8 by /?sb- In Sect. 18.21 we applied 
this correction only to known AGNs. We now consider the pos- 
sibility that galaxies with neither X-ray nor optical evidence of 
an AGN may harbor a dust-obscured AGN whose signature is 
masked by the co-existence of a starburst. 

In the upper panel of Fig. |25]we present the IRS ratios of 
GOODS-Herschel galaxies as in the bottom panel of the right- 
hand part of Fig. |6] The central value and width of the Gaussian 
distribution of main sequence galaxies (Fig. |9] and Eq. |5]i is 
shown with plain and dashed lines respectively. Only 2 % of 
the galaxies fall below the lower limit of the IR main sequence. 
After correcting /7?8 by Rsb, as in Sect. 18.21 we find that the 
fraction of galaxies falling below the IRS-2.4 lower limit of 
the main sequence is increased by a factor of 8, reaching 17 %, 
whereas the number of sources above the main sequence was re- 
duced by 15 %. The effect is stronger in (U)LIRGs, i.e., above 
L8~3xlO'° Lq, where the fraction of galaxies below //?8=2.4 
reaches 25 %. These galaxies present a starburstiness coefficient 
that would normally put them in the high-/7?8 tail of the distri- 
bution. Instead, they exhibit normal values of //?8, and fall down 
to low values of IRE/Rsb- This suggests that part of their 8 pm 
rest-frame radiation is powered by a dust-obscured AGN that 



was not identified from X-ray or optical signatures. These can- 
didate obscured AGN behave similarly to power-law AGN, but 
they were not identified as such because of the presence of a star- 
burst. They are very good candidates for the missing Compton 
Thick AGN needed to explain the peak emission of the cosmic 
X-ray background around 30 keV (Gilli, Comastri & Hasinger 
2007). 

Finally, we applied this technique to a sample of sub- 
millimeter galaxies (SMGs) from Menendez-Delmestre et al. 
(2009) and Pope et al. (2008a). The IRd, and Rsb values of 
these galaxies are shown in Fig. [19] While most objects fol- 
low the trend defined for local galaxies (albeit with a wide dis- 
persion), 1 1 out of a total of 28 SMGs exhibit a starburstiness 
(Rsb) higher than expected for their IRS. Hence if we correct 
IRS for starburstiness in these galaxies, we find evidence for 
the presence of hidden AGN activity. The most extreme cases 
are the SMGs known as "CI", "GN39a" and "GN39b" from 
Pope et al. (2008a). CI, or SMM J123600+621047,is az~2.002 
SMG which has the strongest mid-IR continuum and weakest 
PAH emission lines in a sample of SMGs with Spitzer IRS spec- 
troscopy analyzed by Pope et al. (2008a). Those authors interpret 
this as evidence that 80 % of the mid-IR emission from this ob- 
ject arises from an AGN. It is undetected in the 2 Ms CDF-N 
data (Alexander et al. 2003), and its X-ray to 6 pm luminosity 
ratio indicates that it hosts a Compton-thick AGN (Alexander 
et al. 2008). GN39a and GN39b are both classified as obscured 
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Fig. 25. IR^=L}^/Lfi,ft„ ratio as a function of Lg^m for star- 
forming galaxies, excluding recognized X-ray/optical AGN and 
infrared power-law AGN. The plain and dashed horizontal lines 
show the center and width of the Gaussian distribution of main 
sequence galaxies (Fig. |9] and Eq.|5]l. Upper panel: position of 
the star-forming galaxies as measured. Most fall on the same 
sequence as local star-forming galaxies. Bottom panel: posi- 
tion of the same galaxies after correcting their //?8 ratio for 
starburstiness (unless 0.5</?sb^2, where 7?sb is consistent with 
unity within the error bars). 



AGN based on their strong X-ray hardness ratios, although their 
mid-IR spectra only show 10-40 % AGN contribution. Hence the 
technique appears to be efficient even in the case of extreme sys- 
tems such as distant SMGs. 



9. Discussion and conclusion 

The mode in which galaxies form their stars seems to follow 
some fairly simple scaling laws. The Schmidt-Kennicutt law, 
which connects the surface densities of gas and star formation 
in the local Universe (Kennicutt 1998b), has been recently ex- 
tended to the study of distant galaxies. Two star formation modes 
have thus been identified: so-called "normal" star formation, and 
an accelerated mode, where the star formation efficiency (SFE) 
is increased, probably due to the merger of two galaxies (Daddi 
et al. 2010, Genzel et al. 2010). The different SFEs of these two 
modes are difficult to recognize because of the observational 
challenges associated with measuring the mass and density of 
molecular gas at high redshift. Indeed, the CO luminosity to H2 
conversion factor is poorly known and is based on '^CO 7=1- 
emission locally, while observations of distant galaxies rely 
on higher-7 transitions (see, e.g., Ivison et al. 2011). Similarly, 



star-forming galaxies follow another scaling law: the SFR - M, 
relation, which measures the characteristic time to double the 
stellar mass of a galaxy. At each cosmic epoch, one can iden- 
tify a typical sSFR for star-forming galaxies. This distinguishes 
a main sequence of "normal" star-forming galaxies from a mi- 
nority population of starburst galaxies with elevated sSFR. 

Our analysis of the deep surveys carried out in the open time 
key program GOODS-Herschel allowed us to establish a third 
scaling law for star-forming galaxies relating the total IR lumi- 
nosity of galaxies, Lir, hence their SFR, to the broadband 8 /im 
luminosity, Lg- We showed that the 8 fim bolometric correction 
factor, IRE=LiR/Lii, exhibits a Gaussian distribution containing 
the vast majority of star-forming galaxies both locally and up 
to z~2.5, centered on IR^~4. This defines an IR main sequence 
for star-forming galaxies. Outliers from this main sequence pro- 
duce a tail skewed toward higher values of IRE. We find that this 
sub-population (<20 %) is due to galaxies experiencing compact 
star-formation in a starburst mode. 

The projected star-formation densities of present-day galax- 
ies were estimated from their IR surface brightnesses, Eir, as 
measured from the size of their radio and/or 13.2yum contin- 
uum emission. For distant galaxies, we stacked rest-frame UV- 
2700 A images from HST-ACS in the B, V and / filters for galax- 
ies located at z~ 0.7, 1.2 and 1.8 respectively. We find that at 
all times the projected star-formation density of galaxies in the 
high-IRS tail is more compact (2ir>3x10'° Lokpc"^ at z~0) 
than in galaxies belonging to the IR main sequence, which in- 
cludes distant (U)LIRGs as well. 

Using the more accurate SFRs derived from Herschel data 
for galaxies at 0<z<3, we established the evolution of the typi- 
cal specific SFR (sSFR=SFR/M») for star-forming galaxies. This 
allowed us to separate main sequence and starbursting galaxies 
thanks to a new parameter, labeled "starburstiness" (Rsb), which 
measures the excess sSFR with respect to the SFR - M, main 
sequence, i.e., Rsb=sSFR/sSFRms(z)- We find that galaxies be- 
longing to this main sequence (Rsb = 1±2) also belong to the one 
defined by the Gaussian distribution of IR%, and that the com- 
pact, star-forming galaxies that make up the high-//?8 tail fall 
systematically above the SFR - relation, with strong star- 
burstiness (RsB>2-3). Indeed, we find that IRE is strongly corre- 
lated with RsB in general. Hence IRE appears to be a good proxy 
for identifying compact starbursts, most probably triggered by 
merger events. In the present-day Universe, most (U)LIRGs are 
found to be experiencing compact star-formation during a star- 
burst phase, which is not the case for most distant (U)LIRGs. 
Most probably, the very high SFRs of local (U)LIRGs can only 
be achieved during mergers, whereas distant galaxies are more 
gas-rich and can sustain these large SFRs in other ways. As a 
result of this difference, previous studies that have used local 
(U)LIRG SED templates, with their large, starbursting IRE ra- 
tios, to extrapolate from MIPS 24 yum photometry of galaxies at 
z > 1 .5 have overestimated their total infrared luminosities and 
SFRs, resulting in the so-called "mid-IR excess" issue. 

Using ^-correction as a spectrophotometric tool for convert- 
ing broadband photometric measurements at various redshifts 
into a medium resolution IR SED, we were able to determine 
the prototypical IR SED of MS and SB galaxies with a resolu- 
tion of /l/A/l=25 and 10 respectively. The SED of MS galaxies 
presents strong PAH emission line features, a broad far-IR bump 
resulting from a combination of emission from dust at different 
temperatures ranging typically from 15 to 50 K, and an effective 
dust temperature of 3 1 K, as derived from the peak wavelength 
of the IR SED. Galaxies that inhabit the SB regime instead ex- 
hibit weak PAH equivalent widths and a sharper far-IR bump 
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with an effective dust temperature of ~40 K. Although PAHs are 
stronger in MS than in SB galaxies, they are found to be present 
in the SEDs of both, implying that the IR emission in both pop- 
ulations is primarily powered by star formation and not AGN 
activity. 

Finally, we present evidence that the mid-to-far IR emission 
of X-ray active galactic nuclei is dominantly produced by star 
formation, and that power-law AGNs systematically occur in 
compact, dusty starbursts. After correcting for the excess //?8 
due to star formation - estimated from the starburstiness Rsb 
which we showed correlates with //?8 - we identify candidate 
members of for a sub-population of extremely obscured AGN 
that have not been identified as such by any other method. 

Future studies will be dedicated to understanding the origin 
of the increase of IRE with compactness and starburstiness, by 
separating the relative contributions of PAH lines and contin- 
uum emission, relating the starburstiness with the local environ- 
ment of galaxies as well as with their dust temperature. ALMA 
and eMERLIN will soon provide powerful tools to measure the 
spatial distribution of star formation in distant galaxies at high 
angular resolution, making it possible not only to understand the 
compactness but also the dumpiness of the star-forming regions. 

Spectro-imaging with integral field instruments like 
VLT/SINFONI will also be essential for measuring kinematic 
signatures to assess whether galaxy interactions play an impor- 
tant role in this process (see, e.g., Forster Schreiber et al. 2009, 
Shapiro et al. 2009). Most star formation takes place among 
main sequence galaxies, which suggests that internal, secular 
processes dominate over the role of mergers. This also accounts 
more naturally for the long duty cycle of their star-forming 
phase (Noeske et al. 2007, Daddi et al. 2007a, 2010). Finally, 
we note that this technique of separating MS and SB galaxies 
based on their IR star formation compactness will be extremely 
useful to extrapolate accurate total IR luminosities and SFRs of 
distant galaxies when using the mid-IR camera MIRI on board 
the James Webb Space Telescope. 
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